UNIFORM RANDOM SAMPLING OF PLANAR GRAPHS 
IN LINEAR TIME 



ERIC FUSY 

Abstract. This article introduces new algorithms for the uniform random gen- 
eration of labelled planar graphs. Its principles rely on Boltzmann samplers, as 
recently developed by Duchon, Flajolet, Louchard, and Schaeffer. It combines 
the Boltzmann framework, a suitable use of rejection, a new combinatorial bi- 
jection found by Fusy, Poulalhon and Schaeffer, as well as a precise analytic 
description of the generating functions counting planar graphs, which was re- 
cently obtained by Gimenez and Noy. This gives rise to an extremely efficient 
algorithm for the random generation of planar graphs. There is a preprocessing 
step of some fixed small cost; and the expected time complexity of generation 
is quadratic for exact-size uniform sampling and linear for approximate-size 
sampling. This greatly improves on the best previously known time complexity 
for exact-size uniform sampling of planar graphs with n vertices, which was a 
little over 0{n7). 

This is the extended and revised journal version of a conference paper with 
the title "Quadratic exact-size and linear approximate-size random generation 
of planar graphs", which appeared in the Proceedings of the International Con- 
ference on Analysis of Algorithms (AofA'05), 6-10 June 2005, Barcelona. 



1. Introduction 

A graph is said to be planar if it can be embedded in the plane so that no two edges 
cross each other. In this article, we consider planar graphs that are labelled, i.e., the n 
vertices bear distinct labels in and simple, i.e., with no loop nor multiple edges. 

Statistical properties of planar graphs have been intensively studied [SI HH HO]- Very 
recently, Gimenez and Noy [20] have solved exactly the difficult problem of the asymptotic 
enumeration of labelled planar graphs. They also provide exact analytic expressions for 
the asymptotic probability distribution of parameters such as the number of edges and 
the number of connected components. However many other statistics on random planar 
graphs remain analytically and combinatorially intractable. Thus, it is an important 
issue to design efficient random samplers in order to observe the (asymptotic) behaviour 
of such parameters on random planar graphs. Moreover, random generation is useful 
to test the correctness and efficiency of algorithms on planar graphs, such as planarity 
testing, embedding algorithms, procedures for finding geometric cuts, and so on. 

Denise, Vasconcellos, and Welsh have proposed a first algorithm for the random gen- 
eration of planar graphs [5| , by defining a Markov chain on the set Qn of labelled planar 
graphs with n vertices. At each step, two different vertices v and v' are chosen at ran- 
dom. If they are adjacent, the edge [v, v') is deleted. If they are not adjacent and if the 
operation of adding {v,v') does not break planarity, then the edge {v,v') is added. By 
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Figure 1. Complexities of the random samplers of planar graphs (O* 
stands for a big O taken up to logarithmic factors). 



symmetry of the transition matrix of the Markov chain, the probability distribution con- 
verges to the uniform distribution on Qn- This algorithm is very easy to describe but more 
difficult to implement, as there exists no simple linear-time planarity testing algorithm. 
More importantly, the rate of convergence to the uniform distribution is unknown. 

A second approach for uniform random generation is the recursive method introduced 
by Nijenhuis and Wilf [53] and formalised by Flajolet, Van Cutsem and Zimmermann [15) . 
The recursive method is a general framework for the random generation of combinatorial 
classes admitting a recursive decomposition. For such classes, producing an object of the 
class uniformly at random boils down to producing the decomposition tree corresponding 
to its recursive decomposition. Then, the branching probabilities that produce the de- 
composition tree with suitable (uniform) probability are computed using the coefficients 
counting the objects involved in the decomposition. As a consequence, this method re- 
quires a preprocessing step where large tables of large coefficients are calculated using the 
recursive relations they satisfy. 

Bodirsky et al. have described in [5] the first polynomial-time random sampler for 
planar graphs. Their idea is to apply the recursive method of sampling to a well known 
combinatorial decomposition of planar graphs according to successive levels of connectiv- 
ity, which has been formalised by Tutte |33j . Precisely, the decomposition yields some 
recurrences satisfied by the coefficients counting planar graphs as well as subfamilies (con- 
nected, 2-connected, 3-connected) , which in turn yield an explicit recursive way to generate 
planar graphs uniformly at random. As the recurrences are rather involved, the complex- 
ity of the preprocessing step is large. Precisely, in order to draw planar graphs with n 
vertices (and possibly also a fixed number m of edges) , the random generator described 
in [S] requires a preprocessing time of order O (n^(logn)^(loglogn)) and an auxiliary 
memory of size 0{n^ logn). Once the tables have been computed, the complexity of each 
generation is 0(n'^). A more recent optimisation of the recursive method by Denise and 
Zimmermann [5] — based on controlled real arithmetics — should be applicable; it would 
improve the time complexity somewhat, but the storage complexity would still be large. 

In this article, we introduce a new random generator for labelled planar graphs, which 
relies on the same decomposition of planar graphs as the algorithm of Bodirsky et al. The 
main difference is that wc translate this decomposition into a random generator using the 
framework of Boltzmann samplers, instead of the recursive method. Boltzmann samplers 
have been recently developed by Duchon, Flajolet, Louchard, and Schacffer in [TT] as a 
powerful framework for the random generation of decomposable combinatorial structures. 
The idea of Boltzmann sampling is to gain efficiency by relaxing the constraint of exact-size 
sampling. As we will see, the gain is particularly significant in the case of planar graphs, 
where the decomposition is more involved than for classical classes, such as trees. Given 
a combinatorial class, a Boltzmann sampler draws an object of size n with probability 
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proportional to (or proportional to a;"/n! for labelled objects), where a; is a certain 
real parameter that can be appropriately tuned. Accordingly, the probability distribution 
is spread over all the objects of the class, with the property that objects of the same 
size have the same probability of occurring. In particular, the probability distribution is 
uniform when restricted to a fixed size. Like the recursive method, Boltzmann samplers 
can be designed for any combinatorial class admitting a recursive decomposition, as there 
are explicit sampling rules associated with each classical construction (Sum, Product, Set, 
Substitution). The branching probabilities used to produce the decomposition tree of a 
random object are not based on the coefficients as in the recursive method, but on the 
values at x of the generating functions of the classes intervening in the decomposition. 

In this article, we translate the decomposition of planar graphs into Boltzmann samplers 
and obtain very efhcient random generators that produce planar graphs with a fixed 
number of vertices or with fixed numbers of vertices and edges uniformly at random. 
Furthermore, our samplers have an approximate-size version where a small tolerance, say 
a few percents, is allowed for the size of the output. For practical purpose, approximate- 
size random sampling often suffices. The approximate-size samplers we propose arc very 
efficient as they have linear time complexity. 

Theorem 1 (Samplers with respect to number of vertices). Let n £ IS! be a target size. 
An exact-size sampler 2l„ can be designed so as to generate labelled planar graphs with n 
vertices uniformly at random. For any tolerance ratio e > 0, an approximate-size sampler 
2l„.e can be designed so as to generate planar graphs with their number of vertices in [n{l — 
e), n{l + e)], and following the uniform distribution for each size k € [n{l — e), n(l -t- e)]. 

Under a real- arithmetics complexity model, Algorithm 2l„ is of expected complexity 
0{n^), and Algorithm 2l„.e is of expected complexity 0{n/e). 

Theorem 2 (Samplers with respect to the numbers of vertices and edges). Let n £ N be a 

target size and fi € (1, 3) be a parameter describing the ratio edges-vertices. An exact-size 
sampler 2t„_^ can be designed so as to generate planar graphs with n vertices and [/^nj 
edges uniformly at random. For any tolerance-ratio e > 0, an approximate-size sampler 
^n,ii,e can be designed so as to generate planar graphs with their number of vertices in 
[n(l — e), n(l -|- e)] and their ratio edges/vertices in [fj,(l — e), fi{l -f- e)], and following the 
uniform distribution for each fixed pair (number of vertices, number of edges). 

Under a real- arithmetics complexity model, for a fixed n G (1,3), Algorithm is of 
expected complexity 0^{n^^'^). For fixed constants fi G (1,3) and e > 0, Algorithm 2ln,/j,e 
is of expected complexity Ofj_{n/e) (the bounding constants depend on fi). 

The samplers are completely described in Section 16.11 and Section 16.21 The expected 
complexities will be proved in Section [8l For the sake of simplicity, we give big O bounds 
that might depend on /i and we do not care about quantifying the constant in the big O 
in a precise way. However we strongly believe that a more careful analysis would allow 
us to have a uniform bounding constant (over fi G (1,3)) of reasonable magnitude. This 
means that not only the theoretical complexity is good but also the practical one. (As 
we review in Section [7l we have implemented the algorithm, which easily draws graphs of 
sizes in the range of 10^.) 

Complexity model. Let us comment on the model we adopt to state the complexities of 
the random samplers. We assume here that we are given an oracle, which provides at unit 
cost the exact evaluations of the generating functions intervening in the decomposition 
of planar graphs. (For planar graphs, these generating functions are those of families 
of planar graphs of different connectivity degrees and pointed in different ways.) This 
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assumption, called the "oracle assumption" , is by now classical to analyse the complexity 
of Boltzmann samplers, see for a more detailed discussion; it allows us to separate the 
combinatorial complexity of the samplers from the complexity of evaluating the generating 
functions, which resorts to computer algebra and is a research project on its own. Once the 
oracle assumption is done, the scenario of generation of a Boltzmann sampler is typically 
similar to a branching process; the generation follows a sequence of random choices — 
typically coin flips biased by some generating function values — that determine the shape 
of the object to be drawn. According to these choices, the object (in this article, a planar 
graph) is built effectively by a sequence of primitive operations such as vertex creation, 
edge creation, merging two graphs at a common vertex... The combinatorial complexity 
is precisely defined as the sum of the number of coin flips and the number of primitive 
operations performed to build the object. The (combinatorial) complexity of our algorithm 
is compared to the complexities of the two preceding random samplers in Figure [T] 

Let us now comment on the preprocessing complexity. The implementation of 2l„^e and 
2t„, as well as ^n,p,,e and 2l„^^, requires the storage of a fixed number of real constants, 
which are special values of generating functions. The generating functions to be evaluated 
are those of several families of planar graphs (connected, 2-connected, 3-connected). A 
crucial result, recently established by Gimenez and Noy [20j, is that there exist exact ana- 
lytic equations satisfied by these generating functions. Hence, their numerical evaluation 
can be performed efficiently with the help of a computer algebra system; the complexity 
we have observed in practice (doing the computations with Maple) is of low polynomial 
degree k in the number of digits that need to be computed. (However, there is not yet 
a complete rigorous proof of the fact, as the Boltzmann parameter has to approach the 
singularity in order to draw planar graphs of large size.) To draw objects of size n, the 
precision needed to make the probability of failure small is typically of order log(??) dig- 
itfl Thus the preprocessing step to evaluate the generating functions with a precision of 
log(7i) digits has a complexity of order log(7i)'' (again, this is yet to be proved rigorously). 
The following informal statement summarizes the discussion; making a theorem of it is 
the subject of ongoing research (see the recent article [26]): 

Fact. With high probability, the auxiliary memory necessary to generate planar graphs of 
size n is of order 0(log(n)) and the preprocessing time complexity is of order 0(log(n)'^) 
for some low integer k. 

Implementation and experimental results. We have completely implemented the ran- 
dom samplers stated in Theorem [T] and Theorem [2] Details are given in Section [71 as well 
as experimental results. Precisely, the evaluations of the generating functions of planar 
graphs have been carried out with the computer algebra system Maple, based on the an- 
alytic expressions given by Gimenez and Noy [20] . Then, the random generator has been 
implemented in Java, with a precision of 64 bits for the values of generating functions 
("double" type). Using the approximate-size sampler, planar graphs with size of order 
100,000 are generated in a few seconds with a machine clocked at IGHz. In contrast, the 
recursive method of Bodirsky et al is currently limited to sizes of about 100. 

Having the random generator implemented, we have performed some simulations in or- 
der to observe typical properties of random planar graphs. In particular we have observed 
a sharp concentration for the proportion of vertices of a given degree k in a, random planar 
graph of large size. 



Notice that it is possible to achieve perfect uniformity by caUing adaptive precision routines in case 
of failure, see Denise and Zimmermann [9] for a detailed discussion on similar problems. 
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2. Overview 



The algorithm we describe rches mainly on two ingredients. The first one is a recent 
correspondence, called the closure-mapping, between binary trees and (edge-rooted) 3- 
connected planar graphs |18| , which makes it possible to obtain a Boltzmann sampler for 
3-connected planar graphs. The second one is a decomposition formalised by Tutte [33], 
which ensures that any planar graph can be decomposed into 3-connected components, 
via connected and 2-connccted components. Taking advantage of Tutte's decomposition, 
we explain in Section 3] how to specify a Boltzmann sampler for planar graphs, denoted 
TQ{x, y), from the Boltzmann sampler for 3-connected planar graphs. To do this, we have 
to extend the collection of constructions for Boltzmann samplers, as detailed in and 
develop new rejection techniques so as to suitably handle the rooting/unrooting operations 
that appear alongside Tutte's decomposition. 

Even if the Boltzmann sampler TQ{x, y) already yields a polynomial-time uniform ran- 
dom sampler for planar graphs, the expected time complexity to generate a graph of size 
n (n vertices) is not good, due to the fact that the size distribution of TQ{x^y) is too 
concentrated on objects of small size. To improve the size distribution, we point the ob- 
jects, in a way inspired by [llj . which corresponds to a derivation (differentiation) of the 
associated generating function. The precise singularity analysis of the generating func- 
tions of planar graphs, which has been recently done in |20j . indicates that we have to 
take the second derivative of planar graphs in order to get a good size distribution. In 
Section m we explain how the derivation operator can be injected in the decomposition of 
planar graphs. This yields a Boltzmann sampler TQ" [x, y) for "bi-derived" planar graphs. 
Our random generators for planar graphs are finally obtained as targetted samplers, which 
call TQ"{x,y) (with suitably tuned values of x and y) until the generated graph has the 
desired size. The time complexity of the targetted samplers is analysed in Section [8l This 
eventually yields the complexity results stated in Theorems [T] and [2l The general scheme 
of the planar graph generator is shown in Figure [2l 



In this section, we define Boltzmann samplers and describe the main properties which 
we will need to handle planar graphs. In particular, we have to extend the framework to 
the case of mixed classes, meaning that the objects have two types of atoms. Indeed the 
decomposition of planar graphs involves both (labelled) vertices and (unlabelled) edges. 
The constructions needed to formulate the decomposition of planar graphs are classical 
ones in combinatorics: Sum, Product, Set, Substitutions [3l[14]. In Section [3T2l for each of 
the constructions, we describe a sampling rule, so that Boltzmann samplers can be assem- 
bled for any class that admits a decomposition in terms of these constructions. Moreover, 
the decomposition of planar graphs involves rooting/unrooting operations, which makes 
it necessary to develop new rejection techniques, as described in Section [3.4.31 

3.1. Definitions. A combinatorial class C is a family of labelled objects (structures), 
that is, each object is made of n atoms that bear distinct labels in In addition, the 

number of objects in any fixed size n is finite; and any structure obtained by relabelling 
a structure in C is also in C. The exponential generating function of C is defined as 



3. Boltzmann samplers 
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Figure 2 . The chain of constructions from binary trees to planar graphs. 



where I7I is the size of an object 7 G C (e.g., the number of vertices of a graph). The 
radius of convergence of C (x) is denoted by p. A positive value x is called admissible if 
X G (0, p) (hence the sum defining C{x) converges if x is admissible). 

Boltzmann samplers, as introduced and developed by Duchon et al. in constitute 
a general and efHcient framework to produce a random generator for any decomposable 
combinatorial class C. Instead of fixing a particular size for the random generation, objects 
arc drawn under a probability distribution spread over the whole class. Precisely, given 
an admissible value for C{x), the Boltzmann distribution assigns to each object of C a 
weight 



Notice that the distribution is uniform, i.e., two objects with the same size have the same 
probability to be chosen. A Boltzmann sampler for the labelled class C is a procedure 
TC{x) that, for each fixed admissible x, draws objects of C at random under the distri- 
bution Px- The authors of [TT] give sampling rules associated to classical combinatorial 
constructions, such as Sum, Product, and Set. (For the unlabelled setting, we refer to the 
more recent article [12], and to [4] for the specific case of plane partitions.) 

In order to translate the combinatorial decomposition of planar graphs into a Boltz- 
mann sampler, we need to extend the framework of Boltzmann samplers to the bivariate 
case of mixed combinatorial classes. A mixed class C is a labelled combinatorial class 
where one takes into account a second type of atoms, which are unlabelled. Precisely, an 
object in C = Un,mCn,m has n "labelled atoms" and to "unlabelled atoms", e.g., a graph 
has n labelled vertices and m unlabelled edges. The labelled atoms are shortly called 
L-atoms, and the unlabelled atoms are shortly called U-atoms. For 7 S C, we write I7I for 
the number of L-atoms of 7, called the L-size of 7, and ||7|| for the number of U-atoms 
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of 7, called the U-size of 7. The associated generating function C{x,y) is defined as 

For a fixed real value y > 0, we denote by pciu) the radius of convergence of the function 
X I— > C{x,y). A pair {x,y) is said to be admissible \i x £ (0, pc{y)), which implies that 
'^jec f^lTy"''^" converges and that C{x, y) is well defined. Given an admissible pair {x, y), 
the mixed Boltzmann distribution is the probability distribution Px,j, assigning to each 
object 7 £ C the probability 

P (ry\^ \ ^„ll7ll 

^ x,y\ 1 — \ I II W 

C{x,y) |7|! 

An important property of this distribution is that two objects with the same size-parameters 
have the same probability of occurring. A mixed Boltzmann sampler at (x, y) — shortly 
called Boltzmann sampler hereafter — is a procedure TC{x, y) that draws objects of C at 
random under the distribution ^x,v Notice that the specialization ?/ = 1 yields a classical 
Boltzmann sampler for C. 



3.2. Basic classes and constructions. We describe here a collection of basic classes 
and constructions that are used thereafter to formulate a decomposition for the family of 
planar graphs. 

The basic classes we consider are: 

• The 1-class, made of a unique object of size (both the L-size and the U-size are 
equal to 0), called the 0-atom. The corresponding mixed generating function is 
C{x,y) = l. 

• The L-unit class, made of a unique object that is an L-atom; the corresponding 
mixed generating function is C{x,y) = x. 

• The U-unit class, made of a unique object that is a U-atom; the corresponding 
mixed generating function is C{x, y) = y. 

Let us now describe the five constructions that are used to decompose planar graphs. 
In particular, we need two specific substitution constructions, one at labelled atoms that 
is called L-substitution, the other at unlabelled atoms that is called U-substitution. 

Sum. The sum C := A + B of two classes is meant as a disjoint union, i.e., it is the union 
of two distinct copies of A and B. The generating function of C satisfies 

C{x,y) = A{x,y)+B{x,y). 



Product. The partitional product of two classes A and B is the class C := A-kB of objects 
that are obtained by taking a pair 7 = (71 € 72 € B) and relabelling the L-atoms so 
that 7 bears distinct labels in [I..I7I]. The generating function of C satisfies 

C{x,y) = A{x,y) ■ B{x,y). 



Set>d. For d> and a class B having no object of size 0, any object in C := SET>d{B) is 
a finite set of at least d objects of B, relabelled so that the atoms of 7 bear distinct labels 
in [1 . . I7I]. For d — 0, this corresponds to the classical construction Set. The generating 
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function of C satisfies 

C{x,y) = exp^j^{B{x,y)), where exp>rf(z) := ^ ^. 

k>d 

L-substitution. Given A and B two classes such that B has no object of size 0, the class 
C = ^ S is the class of objects that are obtained as follows: take an object p G A 
called the core-object, substitute each L-atom w of p by an object 7.1, G B, and relabel 
the L-atoms of U^7i, with distinct labels from 1 to |7„|. The generating function of C 
satisfies 

Cix.y) = A{Bix,y),y). 

U-substitution. Given A and B two classes such that B has no object of size 0, the 
class C ^ Aojj B is the class of objects that are obtained as follows: take an object p £ A 
called the core-object, substitute each U-atom e of p by an object 7e € B, and relabel the 
L-atoms of p U (Ue7e) with distinct labels from 1 to \p\ + J^e It^I- assume here that 
the U-atoms of an object of A arc distinguishable. In particular, this property is satisfied 
if ^ is a family of labelled graphs with no multiple edges, since two different edges are 
distinguished by the labels of their extremities. The generating function of C satisfies 

C{x,y) = A{x,B{x,y)). 



3.3. Sampling rules. A nice feature of Boltzmann samplers is that the basic combi- 
natorial constructions (Sum, Product, Set) give rise to simple rules for assembling the 
associated Boltzmann samplers. To describe these rules, we assume that the exact values 
of the generating functions at a given admissible pair {x, y) are known. We will also need 
two well-known probability distributions. 

A random variable follows a Bernoulli law of parameter p e (0, 1) if it is equal to 
1 (or true) with probability p and equal to (or false) with probability 1 — p. 
Given A G IR.+ and d G the conditioned Poisson law Pois>d(A) is the proba- 
bility distribution on 7i>d defined as follows: 

nk) = ^5^;^^' where exp>,(z) :=. . 

For d = 0, this corresponds to the classical Poisson law, abbreviated as Pois. 

Starting from combinatorial classes A and B endowed with Boltzmann samplers r^(a;, y) 
and VB{x, y), Figure [3] describes how to assemble a sampler for a class C obtained from A 
and B (or from A alone for the construction SET>(i) using the five constructions described 
in this section. 

Proposition 3. Let A and B be two mixed combinatorial classes endowed with Boltzmann 
samplers TA{x,y) andTB{x,y). For each of the five constructions {-(-, SET>d, L-subs, 
U-subs}, the sampler TC{x,y), as specified in Figure\3i is a valid Boltzmann sampler for 
the combinatorial class C. 

Proof. 1) Sum: C = A + B. An object of A has probability a(z y) f^^^^^^ definition 
of TA{x,y)) multiplied by (because of the Bernoulli choice) of being drawn by 

TC{x,y). Hence, it has probability y) J^y^^'^^^ '^^ being drawn. Similarly, an object 
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return the L-atom 


U-unit class -Zc/ 


return the U-atom 
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Boltzmann sampler 


Sum e = A + B 


rC{x, y): if Bern (^) , return rA{x, y) 
else return VB(x, y) 


Product C = A*B 


VC{x,y): 7 <— {TA{x,y),VB{x,y)) {independent calls} 
Distributee ABELS (7); return 7 


Set>d C = SET>d(B) 


VC{x,y): k^Vois>d{B{x,y)) 

7 <— {rB{x, y), . . . , VB{x, y)) [k ind. calls} 
DistributeLabels(7); return 7 


L-subs C = AolB 


TC{x,y): 'y^rA{B{x,y),y) 

for each L-atom i; 6 7 do 

replace u by 7„ <— TB{x,y) od {ind. calls} 
DistributeLabels(7); return 7 


U-subs C = AouB 


TC(x,yy. 'y ^rAix,B{x,y)) 

for each U-atom e e 7 do 

replace e by 7e ^ rB{x,y) od {ind. calls} 
DistributeLabels(7); return 7 



Figure 3. The sampling rules associated with the basic classes and the 
constructions. For each rule involving partitional products, there is a re- 
labelling step performed by an auxiliary procedure DistributeLabels. 
Given an object 7 with its L-atoms ranked from 1 to I7I; DistributeLa- 
BELS(7) draws a permutation a of [I..I7I] uniformly at random and gives 
label a{i) to the atom of rank i. 



of B has probability — tttf?/"''^" of being drawn. Hence TC(x, y) is a valid Boltzmann 
sampler for C. 

2) Product: C = A-kB. Define a generation scenario as a pair (71 € .4,72 € .S), together 
with a function a that assigns to each L-atom in 71 U 72 a label i e [l..|7i| + I72I] in a 
bijective way. By definition, TC{x,y) draws a generation scenario and returns the object 
7 g ^ ★ ;B obtained by keeping the secondary labels (the ones given by DistributeLa- 
bels). Each generation scenario has probability 



of being drawn, the three factors corresponding respectively to TA{x,y), TB{x,y), and 
DistributeLabels(7). Observe that this probability has the more compact form 

1 1 a^l^l ii^ii 

ry ■ 



|7i|!|72|!C(a;,y) |7|!' 

Given 7 G ^ t*: S, let 71 be its first component (in A) and 72 be its second component 
(in B). Any relabelling of the labelled atoms of 71 from 1 to |7i| and of the labelled 
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atoms of 72 from 1 to I72I induces a unique generation scenario producing 7. Indeed, the 
two relabeUings determine unambiguously the relabeUing permutation a of the generation 
scenario. Hence, 7 is produced from |7i|!|72|! different scenarios, each having probabihty 
7i|! | 72|!C(g: v) f^^^^'''^^' ^ consequcncc, 7 is drawn under the Boltzmann distribution. 

3) Set>d- C = SET>rf(S). In the case of the construction SET>d, a generation scenario is 
defined as a sequence (71 € S, . . . , 7fc G S) with k > d, together with a function a that 
assigns to each L-atom in 71 U • • • U 7fe a label i £ [l--|7i| + • • • + |7a;|] in a bijective way. 
Such a generation scenario produces an object 7 G Set>c;(S). By definition of TC{x,y), 
each scenario has probability 



exp>^(i3(a:,2/)) fc! / 1^^^^^ y)|7.|7 (l7i| + • • • + l7fe|)! 



the three factors corresponding respectively to drawing Pois> d(B(x,y)), drawing the se- 
quence, and the relabelling step. This probability has the simpler form 



1 Il7ll 

ry 



fc ^ 

n iTjT- 



k\C{x,y)\jr fil7, 

For k > d, an object 7 e SET>rf(S) can be written as a sequence 71, . . . , 7*; in fc! different 
ways. In addition, by a similar argument as for the Product construction, a sequence 
71, . . . ,7fe is produced from Y['i=i 17*1- different scenarios. As a consequence, 7 is drawn 
under the Boltzmann distribution. 

4) L- substitution: C ~ Aoj^ B. For this construction, a generation scenario is defined as 
a core-object p £ A, a. sequence 71, ... , 7|p| of objects of B {-fi stands for the object of B 
substituted at the atom i of p), together with a function cr that assigns to each L-atom 
in 7i U • • • U 7|p| a label i £ [I..I71 1 -|- • • • -I- |7|p| |] in a bijective way. This corresponds to 
the scenario of generation of an object j E Ao^ B hy the algorithm rC{x,y), and this 
scenario has probability 

i B(x^^||p||^ xMylMl \ I 




AiB{x,y),y) \p\l " J \l\ Bix,y)\j,\l j + ■ ■ ■ + \j\,\\y: 
which has the simpler form 



a;l7lyll7ll 



-n- 

nil li U.-l 



C{x,yM\\p\ll\h,\l- 

Given j G Aoj^ B, labelling the core-object p £ A with distinct labels in and each 

component (7i)i<i<|p| with distinct labels in [l..|7i|] induces a unique generation scenario 
producing 7. As a consequence, 7 is produced from |p|!n'=i 17*1- scenarios, each having 
probability ^^T^ypj^j^ 11 '=1 J^.' Hence, 7 is drawn under the Boltzmann distribution. 

5) U- substitution: C = Aou B. A generation scenario is defined as a core-object p G -4, 
a sequence 71, . . . , 7||p|| of objects of B (upon giving a rank to each unlabelled atom of p, 
7i stands for the object of B substituted at the U-atom of rank z in p), and a function a 
that assigns to each L-atom in p U 71 U • • • U 7||p|| a label i £ [l..|p| + |7i| + ■ • • + l7||p|||]- 
This corresponds to the scenario of generation of an object £ Ao^ B by the algorithm 
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TC{x,y); this scenario has probabihty 

1 xM \ :^l7.|yll7.|| 

A{x,B{x,y))\p\\ ^""'y' ) [l\B{x,y)\^,\\j V(|p| + I71I + • • • + |7||p||I)! 



This expression has the simpler form 



-tt — 

oil 11 u.ir 



C{x,yM\\p\\l\ |7.|! 

Given 7 G Aoir B, labelhng the eore-object p £ A with distinct labels in and 
each component (7i)i<i<||p|| with distinct labels in [l..|7i|] induces a unique generation 
scenario producing 7. As a consequence, 7 is produced from |p|!n!=i 17*1- scenarios, 

T I I I T I I II II 

each having probability y)|7|! 111=1 Hence, 7 is drawn under the Boltzmann 

distribution. □ 



Example. Consider the class C of rooted binary trees, where the (labelled) atoms are the 
inner nodes. The class C has the following decomposition grammar, 

C = (C + 1)*Z*(C + 1). 

Accordingly, the series C{x) counting rooted binary trees satisfies C{x) = x (1 + C(.x)) . 
(Notice that C{x) can be easily evaluated for a fixed real parameter x < pc = 1/4.) 

Using the sampling rules for Sum and Product, we obtain the following Boltzmann 
sampler for binary trees, where {•} stands for a node: 

rC{x) : return (r(l + C)(.t), {•}, r(l + C)(a:)) {independent calls} 
r(l + C){x) : if Bern (jqr^j^) return leaf 
else return TC{x) 

Distinct labels in [l..|7|] might then be distributed uniformly at random on the atoms of 
the resulting tree 7, so as to make it well-labelled (sec Rcmark|4]bclow). Many more exam- 
ples are given in [TT] for labelled (and unlabelled) classes specified using the constructions 
{+,*,Set} □ 

Remark 4. In the sampling rules (Figure [3]), the procedure DistributeLabels(7) 
throws distinct labels uniformly at random on the L-atoms of 7. The fact that the re- 
labelling permutation is always chosen uniformly at random ensures that the process of 
assigning the labels has no memory of the past, hence DistributeLabels needs to be 
called just once, at the end of the generation procedure. (A similar remark is given by 
Flajolet et al. in [T31 Sec. 3] for the recursive method of sampling.) 

In other words, when combining the sampling rules given in Figure [3] in order to de- 
sign a Boltzmann sampler, we can forget about the calls to DistributeLabels, see for 
instance the Boltzmann sampler for binary trees above. In fact, we have included the 
DistributeLabels steps in the definitions of the sampling rules only for the sake of 
writing the correctness proofs (Proposition [3]) in a proper way. 

3.4. Additional techniques for Boltzmann sampling. As the decomposition of pla- 
nar graphs we consider is a bit involved, we need a few techniques in order to properly 
translate this decomposition into a Boltzmann sampler. These techniques, which are 
described in more detail below, are: bijections, pointing, and rejection. 
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3.4.1. Combinatorial isomorphisms. Two mixed classes A and B are said to be isomorphic, 
shortly written a.s A ~ B, if there exists a bijection <I> between A and B that preserves the 
size parameters, i.e., preserves the L-sizc and the U-sizc. (This is equivalent to the fact 
that the mixed generating functions of A and B are equal.) In that case, a Boltzmann 
sampler TA{x, y) for the class A yields a Boltzmann sampler for B via the isomorphism: 
TB{x, y) : 7 ^ r^(a;, y); return $(7). 

3.4.2. L- derivation, U- derivation, and edge-rooting. In order to describe our random sam- 
pler for planar graphs, we will make much use of derivative operators. The L-derived class 
of a mixed class C ~ Un^mCn,m (shortly called the derived class of C) is the mixed class 
C = Un.mCn m of objccts in C whcrc the greatest label is taken out, i.e., the L-atom with 
greatest label is discarded from the set of L-atoms (see the book by Bergeron, Labelle, 
Leroux for more details and examples). The class C can be identified with the pointed 
class C* of C, which is the class of objects of C with a distinguished L-atom. Indeed the 
discarded atom in an object of C plays the role of a pointed vertex. However the impor- 
tant difference between C and C* is that the distinguished L-atom does not count in the 
L-size of an object in C . In other words, C* = * C . Clearly, for any integers n,m, 
Cn_i ,„ identifies to Cn^m, so that the generating function C'{x,y) of C satisfies 

(1) C'ix,y) = IC"-'"!?;^^!)!?^" = d,.Cix,y). 

The U-derived class of C is the class C of objects obtained from objects of C by discarding 
one U-atom from the set of U-atoms; in other words there is a distinguished U-atom that 
does not count in the U-size. As in the definition of the U-substitution, we assume that all 
the U-atoms are distinguishable, for instance the edges of a simple graph are distinguished 
by the labels of their extremities. In that case, |C„ = m\Cn^m\i so that the generating 
function C_{x, y) of C satisfies 

(2) C(x,y) =^m|C„.™|— y™-i =9,C(x,2/). 

For the particular case of planar graphs, we will also consider edge-rooted objects 
(shortly called rooted objects), i.e., planar graphs where an edge is "marked" (distin- 
guished) and directed. In addition, the root edge, shortly called the root, is not counted 
as an unlabelled atom, and the two extremities of the root do not count as labelled atoms 
(i.e., are not labelled). The edge-rooted class of C is denoted by C . Clearly we have 
-k C ~ 2 ^C. Hence, the generating function C {x, y) of C satisfies 

-> 2 

(3) C{x,y) = ^dyC{x,y). 

3.4.3. Rejection. Using rejection techniques offers great fiexibility to design Boltzmann 
samplers, since it makes it possible to adjust the distributions of the samplers. 

Lemma 5 (Rejection). Given a combinatorial class C, letW :C ^ and p : C [0, 1] 
be two functions, called weight-function and rejection-function, respectively. Assume that 
W is summable, i.e., X)76C ^(t) "is finite. Let% be a random generator for C that draws 
each object 7 € C with probability proportional to W{'^). Then, the procedure 

Stroj : repeat 21 ^ 7 until Bern(p(7)); return 7 
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is a random generator on C, which draws each object 7 G C with probability proportional 
to W{"f)p{'y). 

Proof. Define W := J^^ec ^(t)- definition. 21 draws an object 7 e C with probability 
P(7) := W{'^)/W. Let Prej be the probabihty of failure of Slrej at each attempt. The 
probability Prej(7) that 7 is drawn by Strej satisfies Prej(7) = ^'(7)p(7)+PrcjPrej(7), where 
the first (second) term is the probability that 7 is drawn at the first attempt (at a later 
attempt, respectively). Hence, Prej(7) = ^(7)p(7)/(l -Prcj) = W {j)p{'-f) / {W ■ {1 - Prcj)) , 
i.e., Prcj(7) is proportional to ^^(7)^(7). □ 

Rejection techniques are very useful for us to change the way objects are rooted. Typ- 
ically it helps us to obtain a Boltzmann sampler for A' from a Boltzmann sampler for A 
and vice versa. As we will use this trick many times, we formalise it here by giving two 
explicit procedures, one from L-dcrived to U-derivcd objects, the other one from U-derived 
to L-derivcd objects. 



Lderived^Uderived 

INPUT: a mixed class A such that ajj^j^ := sup..yg_4-Lpy is finite, 

a Boltzmann sampler TA'{x, y) for the L-derivcd class A! 

OUTPUT: a Boltzmann sampler for the U-derived class A, defined as: 

r^(a;, y): repeat 7 ^ TA'{x, y) {at this point 7 G A'} 

give label I7I + 1 to the discarded L-atom of 7; 
{so I7I increases by 1, and 7 e yl} 

."[//L I7I 

choose a U-atom uniformly at random and discard it 

from the set of U-atoms; {so ||7|| decreases by 1, and 7 G A\ 
return 7 



until Bern 



Lemma 6. The procedure Lderived^Uderived yields a Boltzmann sampler for the 
class A from a Boltzmann sampler for the class A' . 

Proof. First, observe that the sampler is well defined. Indeed, by definition of the pa- 
rameter ajj/L, the Bernoulli choice is always valid (i.e., its parameter is always in [0, 1]). 
Notice that the sampler 
7^r^'(x,2/); 

give label I7I + 1 to the discarded L-atom of 7; 
return 7 

is a sampler for A that outputs each object j ^ A with probability ^'(^ y) (f^jii)! ^/^^^^^ 
because An,m identifies to A'^^i „. In other words, this sampler draws each object "f £ A 
with probability proportional to |7|^^2/"''"- Hence, according to Lemma[51 the rcpeat- 
until loop of the sampler TA{x, y) yields a sampler for A such that each object has 
probability proportional to II7II j^y"''^"- As each U-atom has probability I/II7II of being 
discarded, the final sampler is such that each object 7 S ^4 has probability proportional 
to j^y"'''"- So TA{x,y) is a Boltzmann sampler for A. □ 



We define a similar procedure to go from a U-derived class to an L-derived class: 
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Uderived— s-LDERIVED 

INPUT: a mixed class A such that a^/^/ := sup^g_^-j-j — rr is finite, 



a Bohzmann sampler TA{x, y) for the U-derived class A 

OUTPUT: a Boltzmamr sampler for the L-derived class A! ^ defined as: 

TA'{x, y): repeat 7 <— TA{x, y) {at this point 7 e A\ 

take the discarded U-atom of 7 back in the set of U-atoms; 
{so II7II increases by 1, and ^ & A} 

mitil Bern | 

\aL/u WlW 

discard the L-atom with greatest label from the set of L-atoms; 
{so I7I decreases by 1, and 76^'} 
return 7 



Lemma 7. The ■procedure Uderived^Lderived yields a Boltzmann sampler for the 
class A' from a Boltzmann sampler for the class A. 

Proof. Similar to the proof of Lemma [51 The sampler TA'{x,y) is well defined, as the 
Bernoulli choice is always valid (i.e., its parameter is always in [0,1]). Notice that the 
sampler 

7 ^ TA{x,y); 

take the discarded U-atom back to the set of U-atoms of 7; 
return 7 

is a sampler for A that outputs each object ^ € A with probability ]4(7^||7l||^y"'''"~^, 
(because an object ^An,m gives rise to m objects in A^ m-i)' i-*^-' with probability propor- 
tional to ||7|||^y"'''"- Hence, according to LcmmaO the rcpeat-until loop of the sampler 
TA! {x, y) yields a sampler for A such that each object 76^ has probability proportional 
to 171^:^2/"'''", i.e., proportional to (p^|3YjT 2/ Hence, by discarding the greatest L-atom 

(i.e., I7I <— I7I — 1), we get a probability proportional to for every object 7 £ A' , 

i.e., a Boltzmann sampler for A! . □ 

Remark 8. We have stated in Remark |4] that, during a generation process, it is more 
convenient in practice to manipulate the shapes of the objects without systematically 
assigning labels to them. However, in the definition of the sampler TAHx^y), one step is 
to remove the greatest label, so it seems we need to look at the labels at that step. In fact, 
as we consider here classes that are stable under relabelling, it is equivalent in practice to 
draw uniformly at random one vertex to play the role of the discarded L-atom. 

4. Decomposition of planar graphs and Boltzmann samplers 

Our algorithm starts with the generation of 3-connected planar graphs, which have 
the nice feature that they are combinatorially tractable. Indeed, according to a theorem 
of Whitney [35], 3-connected planar graphs have a unique embedding (up to reflection), 
so they are equivalent to 3-connected planar maps. Following the general approach in- 
troduced by Schaeffer [59], a bijection has been described by the author, Poulalhon, and 
Schaeffer [18] to enumerate 3-connected maps [18] from binary trees, which yields an 
explicit Boltzmann sampler for (rooted) 3-connected maps, as described in Section [4. II 

The next step is to generate 2-connected planar graphs from 3-connected ones. We 
take advantage of a decomposition of 2-connected planar graphs into 3-connected planar 
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Figure 4. The complete scheme to obtain a Boltzmann sampler for pla- 
nar graphs. The classes are to be defined all along Section ID 



Figure 5. The connected planar graphs with at most four vertices (the 
2-connected ones are surrounded). Below each graph is indicated the 
number of distinct labcUings. 



components, which has been formalised by Trakhtenbrot [31] (and later used by Walsh [34] 
to count 2-connected planar graphs and by Bender, Gao, Wormald to obtain asymptotic 
enumeration [1]). Finally, connected planar graphs are generated from 2-connected ones 
by using the well-known decomposition into blocks, and planar graphs are generated from 
their connected components. Let us mention that the decomposition of planar graphs into 
3-connected components has been completely formalised by Tutte [33j (though we rather 
use here formulations of this decomposition on rooted graphs, as Trakhtenbrot did). 
The complete scheme we follow is illustrated in Figure H] 

Notations. Recall that a graph is fc-connected if the removal of any set of A: — 1 vertices 
does not disconnect the graph. In the sequel, we consider the following classes of planar 
graphs: 

Q: the class of all planar graphs, including the empty graph, 
Qi : the class of connected planar graphs with at least one vertex, 
Q2: the class of 2-connected planar graphs with at least two vertices, 
Q^: the class of 3-conncctcd planar graphs with at least four vertices. 
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All these classes are considered as mixed, with labelled vertices and unlabelled edges, 
i.e., the L-atoms are the vertices and the U-atoms are the edges. Let us give the first 
few terms of their mixed generating functions (see also Figure [51 which displays the first 
connected planar graphs): 

G[x,y) = l + a; + §^(l + y) + f^(l + 3y + 3y2 + y3) + ... 

G2{x,y) = |iy+|iy3 + |l(3y4 + 6y5 + y6) + |i(12y5 + 70y6 + 100y7+15y8 + 10j^9) + ... 

G:i{x,y) = 1^2/*^ + 1^(152/8 + 102/9) + |^(60y9 + 432j/i" + 540y" + 195j/i2) + ... 

Observe that, for a mixed class A of graphs^ the derived class ,4', as defined in Sec- 
tion 13.4.21 is the class of graphs in A that have one vertex discarded from the set of 
L-atoms (this vertex plays the role of a distinguished vertex); A is the class of graph in A 
with one edge discarded from the set of U-atoms (this edge plays the role of a distinguished 
edge); and A is the class of graphs in A with an ordered pair of adjacent vertices (m,w) 
discarded from the set of L-atoms and the edge (u, v) discarded from the set of U-atoms 
(such a graph can be considered as rooted at the directed edge (w, v)). 

4.1. Boltzmann sampler for 3-connected planar graphs. In this section wc develop 
a Boltzmann sampler for 3-conncctcd planar graphs, more precisely for edge-rooted ones, 
i.e., for the class ^3. Our sampler relics on two results. First, we recall the equivalence 
between 3-connected planar graphs and 3-connccted maps, where the terminology of map 
refers to an explicit embedding. Second, we take advantage of a bijection linking the 
families of rooted 3-connected maps and the (very simple) family of binary trees, via 
intermediate objects that are certain quadrangular dissections of the hexagon. Using the 
bijection, a Boltzmann sampler for rooted binary trees is translated into a Boltzmann 
sampler for rooted 3-connected maps. 

4.1.1. Maps. A map on the sphere [planar map, resp.) is a connected planar graph em- 
bedded on the sphere (on the plane, resp.) up to continuous deformation of the surface, 
the embedded graph carrying distinct labels on its vertices (as usual, the labels range from 
1 to n, the number of vertices). A planar map is in fact equivalent to a map on the sphere 
with a distinguished face, which plays the role of the unbounded face. The unbounded 
face of a planar map is called the outer face, and the other faces are called the inner faces. 
The vertices and edges of a planar map arc said to be outer or inner whether they are 
incident to the outer face or not. A map is said to be rooted if the embedded graph is 
edge-rooted. The root vertex is the origin of the root. Classically, rooted planar maps are 
always assumed to have the outer face on the right of the root. With that convention, 
rooted planar maps are equivalent to rooted maps on the sphere (given a rooted map on 
the sphere, take the face on the right of the root as the outer face). See Figure |6l(c) for 
an example of rooted planar map, where the labels are forgotter0. 

4.1.2. Equivalence between 3-connected planar graphs and 3-connected maps. A well known 
result due to Whitney [35] states that a labelled 3-connected planar graph has a unique 
embedding on the sphere up to continuous deformation and reflection (in general a planar 
graph can have many cmbcddings). Notice that any 3-connected map on the sphere with 



Classically, rooted maps arc considered in the literature without labels on the vertices, as the root is 
enough to avoid symmetries. Nevertheless, it is convenient here to keep the framework of mixed classes 
for maps, as we do for graphs. 
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at least 4 vertices differs from its mirror-image, due to the labels on the vertices. Hence 
every 3-connected planar graph with at least 4 vertices gives rise exactly to two maps on 
the sphere. The class of 3-connectcd maps on the sphere with at least 4 vertices is denoted 
by A^3. As usual, the class is mixed, the L-atoms being the vertices and the U-atoms 
being the edges. Whitney's theorem ensures that 

(4) A^3-2*53. 

Here we make use of the formulation of this isomorphism for edge-rooted objects. The 
mixed class of rooted 3-connected planar maps with at least 4 vertices is denoted by 
Ms, where — as for edge- rooted graphs — the L-atoms are the vertices not incident to 
the root-edge and the U-atoms are the edges except the root. Equation (j4]) becomes, for 
edge-rooted objects: 

(5) Ml~2i.gl. 

Thanks to this isomorphism, finding a Boltzmann sampler TQ3{z,w) for edge-rooted 
3-connected planar graphs reduces to finding a Boltzmann sampler rAis{z,w) for rooted 
3-connected maps, upon forgetting the embedding. 



4.1.3. 3-connected maps and irreducible dissections. We consider here some quadrangular 
dissections of the hexagon that are closely related to 3-connected planar maps. (We will 
see that these dissections can be efficiently generated at random, as they are in bijection 
with binary trees.) 

Precisely, a quadrangulated map is a planar map (with no loop nor multiple edges) such 
that all faces except maybe the outer one have degree 4; it is called a quadrangulation 
if the outer face has degree 4. A quadrangulated map is called bicolored if the vertices 
are colored black or white such that any edge connects two vertices of different colors. A 
rooted quadrangulated map (as usual with planar maps, the root has the outer face on 
its right) is always assumed to be endowed with the unique vertex bicoloration such that 
the root vertex is black (such a bicoloration exists, as all inner faces have even degree). 
A quadrangulated map with an outer face of degree more than 4 is called irreducible 
if each 4-cycle is the contour of a face. In particular, we define an irreducible dissec- 
tion of the hexagon — shortly called irreducible dissection hereafter — as an irreducible 
quadrangulated map with an hexagonal outer face, see Figure IHJb) for an example. A 
quadrangulation is called irreducible if it has at least 2 inner vertices and if every 4- 
cycle, except the outer one, delimits a face. Notice that the smallest irreducible dissection 
has one inner edge and no inner vertex (see Figure [7]), whereas the smallest irreducible 
quadrangulation is the embedded cube, which has 4 inner vertices and 5 inner faces. We 
consider irreducible dissections as objects of the mixed type, the L-atoms are the black 
inner vertices and the U-atoms are the inner faces. It proves more convenient to consider 
here the irreducible dissections that are asymmetric, meaning that there is no rotation 
fixing the dissection. The four non-asymmetric irreducible dissections are displayed in 
Figure [7{b), all the other ones are asymmetric either due to an asymmetric shape or due 
to the labels on the black inner vertices. We denote by X the mixed class of asymmetric 
bicolored irreducible dissections. We define also J as the class of asymmetric irreducible 
dissections that carry a root (outer edge directed so as to have a black origin and the 
outer face on its right), where this time the L-atoms are the black vertices except two of 
them (say, the origin of the root and the next black vertex in ccw order around the outer 
face) and the U-atoms are all the faces, including the outer one. Finally, we define Q as 
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the mixed class of rooted irreducible quadrangulations, where the L-atoms are the black 
vertices except those two incident to the outer face, and the U-atoms are the inner faces. 

Irreducible dissections arc closely related to 3-connccted maps, via a classical corre- 
spondence between planar maps and quadrangulations. Given a bicolored rooted quad- 
rangulation k, the primal map of k is the rooted map /i whose vertex set is the set of black 
vertices of k, each face / of k giving rise to an edge of /i connecting the two (opposite) 
black vertices of /, see Figure [Hc)-(d). The map fi is naturally rooted so as to have the 
same root- vertex as ft. 

Theorem 9 (MuUin and Schellenberg |24| ) . The primal-map construction is a bijection 
between rooted irreducible quadrangulations with n black vertices and m faces, and rooted 3- 
connected maps with n vertices and m edge^. In other words, the primal-map construction 
yields the combinatorial isomorphism 

(6) Qc^mI. 

In addition, the construction of a 3-connected map from an irreducible quadrangulation 
takes linear time. 

The link between J' and Ai^ is established via the family Q, which is at the same time 
isomorphic to and closely related to J^. Let k be a rooted irreducible quadrangulation, 
and let e be the edge following the root in cw order around the outer face. Then, deleting e 
yields a rooted irreducible dissection S. In addition it is easily checked that S is asymmetric, 
i.e., the four non-asymmetric irreducible dissections, which are shown in Figure [Tl^b) , can 
not be obtained in this way. Hence the so-called root-deletion mapping is injective from 
Q to j/. The inverse operation — called the root-addition mapping — starts from a rooted 
irreducible dissection 6, and adds an outer edge from the root-vertex of 6 to the opposite 
outer vertex. Notice that the rooted quadrangulation obtained in this way might not be 
irreducible. Precisely, a non-separating 4-cyclc appears iff S has an internal path (i.e., a 
path using at least one inner edge) of length 3 connecting the root vertex to the opposite 
outer vertex. A rooted irreducible dissection S is called admissible iff it has no such path. 
The subclass of rooted irreducible dissections that are admissible is denoted by J7a. We 
obtain the following result, already given in [18| : 

Lemma 10. The root-addition mapping is a bijection between admissible rooted irreducible 
dissections with n black vertices and m faces, and rooted irreducible quadrangulations with 
n black vertices and m inner faces. In other words, the root-addition mapping realises the 
combinatorial isomorphism 

(7) Ja ^ Q. 

To sum up, we have the following link between rooted irreducible dissections and rooted 
3-connected maps: 

Notice that we have a combinatorial isomorphism between J7a and Ms: the root-edge 
addition combined with the primal map construction. For i5 S J7a, the rooted 3-connccted 
map associated with S is denoted Primal(5). 

As we see next, the class I (and also the associated rooted class J') is combinatorially 
tractable, as it is in bijection with the simple class of binary trees; hence irreducible 
dissections are easily generated at random. 



'More generally, the bijection holds between rooted quadrangulations and rooted 2-connectcd maps. 
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(a) (b) (c) (d) 

Figure 6. (a) A binary tree, (b) the associated irreducible dissection S 
(rooted and admissible), (c) the associated rooted irreducible quadran- 
gulation k — Add ((5), (d) the associated rooted 3-connected map /i = 
Primal((5). 




(a) (b) 

Figure 7. (a) The four non-asymmetric bicolored binary trees, (b) The 
four non-asymmetric bicolored irreducible dissections. 



4.1.4. Bijection between binary trees and irreducible dissections. There exist by now sev- 
eral elegant bijections between families of planar maps and families of plane trees that 
satisfy simple context-free decomposition grammars. Such constructions have first been 
described by Schaeffer in his thesis [52] , and many other families of rooted maps have been 
counted in this way [17l [27l [28l [7]. The advantage of bijective constructions over recur- 
sive methods for counting maps [32] is that the bijections yield efficient — linear-time — 
generators for maps, as random sampling of maps is reduced to the much easier task of 
random sampling of trees, see [3D]. The method has been recently applied to the family 
of 3-connected maps, which is of interest here. Precisely, as described in [TH], there is a 
bijection between binary trees and irreducible dissections of the hexagon, which, as we 
have seen, are closely related to 3-connected maps. 

We define an unrooted binary tree, shortly called a binary tree hereafter, as a plane tree 
(i.e., a planar map with a unique face) where the degree of each vertex is either 1 or 3. 
The vertices of degree 1 (3) are called leaves (nodes, resp.). A binary tree is said to be 
bicolored if its nodes are bicolored so that any two adjacent nodes have different colors, 
see FigurelHKa) for an example. In a bicolored binary tree the L-atoms are the black nodes 
and the U-atoms are the leaves. A bicolored binary tree is called asymmetric if there is no 
rotation-symmetry fixing it. Figure [7| displays the four non-asymmetric bicolored binary 
trees; all the other bicolored binary trees are asymmetric, either due to the shape being 
asymmetric, or due to the labels on the black nodes. We denote by K. the mixed class 
of asymmetric bicolored binary trees (the requirement of asymmetry is necessary so that 
the leaves are distinguishable). 

The terminology of binary tree refers to the fact that, upon rooting a binary tree at an 
arbitrary leaf, the neighbours in clockwise order around each node can be classified as a 
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father (the neighbour closest to the root), a right son, and a left son, which corresponds 
to the classical definition of rooted binary trees, as considered in Example 13.31 

Proposition 11 (Fusy, Poulalhon, and Schaeffer [18]). For n > and m > 2, there 
exists an explicit hijection, called the closure-mapping, between bicolored binary trees with 
n black nodes and m leaves, and bicolored irreducible dissections with n black inner nodes 
and m inner faces; moreover the 4 non- asymmetric bicolored binary trees are mapped to 
the 4 non- asymmetric irreducible dissections. In other words, the closure-mapping realises 
the combinatorial isomorphism 

(8) /C ~ I. 

The construction of a dissection from a binary tree takes linear time. 

Let us comment a bit on this bijective construction, which is described in detail in [18j . 
Starting from a binary tree, the closure-mapping builds the dissection face by face, each 
leaf of the tree giving rise to an inner face of the dissection. More precisely, at each 
step, a "leg" (i.e., an edge incident to a leaf) is completed into an edge connecting two 
nodes, so as to "close" a quadrangular face. At the end, an hexagon is created outside of 
the figure, and the leaves attached to the remaining non-completed legs are merged with 
vertices of the hexagon so as to form only quadrangular faces. For instance the dissection 
of Figure (Hub) is obtained by "closing" the tree of Figure [Hlja). 



4.1.5. Boltzmann sampler for rooted bicolored binary trees. We define a rooted bicolored 
binary tree as a binary tree with a marked leaf discarded from the set of U-atoms. Notice 
that the class of rooted bicolored binary trees such that the underlying unrooted binary 
tree is asymmetric is the U-derived class JC. 

In order to write down a decomposition grammar for the class IC — to be translated 
into a Boltzmann sampler — we define some refined classes of rooted bicolored binary trees 
(decomposing ^ is a bit involved since we have to forbid the 4 non- asymmetric binary 
trees): TZ, is the class of black-rooted binary trees (the root leaf is connected to a black 
node) with at least one node, and TZo is the class of white-rooted binary trees (the root leaf 
is connected to a white node) with at least one node. We also define Tzi^^'' (TZ'f^'') as the 
class of black-rooted (white-rooted, resp.) bicolored binary trees such that the underlying 
unrooted binary tree is asymmetric. Hence IC = TZ^'^^ -{-Tzi^^^^ . We introduce two auxiliary 
classes; TZ, is the class of black-rooted binary trees except the (unique) one with one black 
node and two white nodes; and TZo is the class of white-rooted binary trees except the two 
ones resulting from rooting the (unique) bicolored binary tree with one black node and 
three white nodes (the 4th one in Figure [7]J a)), in addition, the rooted bicolored binary 
tree with two leaves (the first one in Figure [7{a)) is also included in the class TZq. 

The decomposition of a bicolored binary tree at the root yields a complete decom- 
position grammar, given in Figure [HI for the class IC — Tzi^^^ -\- TZi^'''^ . This grammar 
translates to a decomposition grammar involving only the basic classes {Zl, Zu} and the 
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7^o = 



n.- 











fas) 

Figure 8. The decomposition grammar for the two classes TZ, and 
Ti!'^'^ of rooted bicolorcd binary trees such that the underlying binary 
tree is asymmetric. 



constructions {+,*} {Zl stands for a black node and Zjj stands for a non-root leaf): 



(9) 



^(as) 



^(as) 



7^, 



(as) 
o 

n. 



TZo * Zl * Zjj + Zjj -k Zl * TZo 

n.-kZu + Zu^n. + nl, 



TZo* Zl* Zfj + Zfj k ZL^lio - 
Zu + Ti,irZu + Zu*Ti,+ Til 
{Zu +TZo)-kZL-kiZu +TZo), 
{Zu +TZ,)*{Zu +TZ,). 



Zl^TZI 



TZo*Zl*TZo 



In turn, this grammar is translated into a Boltzmann sampler TJC(z,w) for the class 
IC using the sampling rules given in Figure [31 similarly as we have done for the (simpler) 
class of complete binary trees in Example 1. 

4.1.6. Boltzmann sampler for bicolored binary trees. We describe in this section a Boltz- 
mann sampler TIC{z,w) for asymmetric bicolored binary trees, which is derived from the 
Boltzmann sampler T!C{x, y) described in the previous section. Observe that each asym- 
metric binary tree in /C„.m gives rise to m rooted binary trees in ^„ each of the 
m leaves, which are distinguishable, might be chosen to be discarded from the set of U- 
atoms. Hence, each object of /C„^,„ has probability K_{z,w)~^mz"- /nly™~^ to be chosen 
when calling TIC{z, w) and taking the distinguished atom back into the set of U-atoms. 
Hence, from the rejection lemma (Lemma [5]), the sampler 
repeat 7 <— T!C{z,w); 

take the distinguished U-atom back into the set of U-atoms; 

{so II7II increases by 1 and now 7 e /C} 

until Bern (y^) ; 

return 7 
is a Boltzmann sampler for IC. 

However, this sampler is not efficient enough, as it uses a massive amount of rejection to 
draw a tree of large size. Instead, we use an early-abort rejection algorithm, which allows 
us to "simulate" the rejection step all along the generation, thus making it possible to 
reject before the entire object is generated. We find it more convenient to use the number 
of nodes, instead of leaves, as the parameter for rejection (the subtle advantage is that 
the generation process TIC{z, w) builds the tree node by node). Notice that the number of 
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leaves in an unrooted binary tree 7 is equal to 2 + N^j), with A^(7) the number of nodes 
of 7. Hence, the rejection step in the sampler above can be replaced by a Bernoulli choice 
with parameter 2/{N{'y) + 2). We now give the early-abort algorithm, which repeats 
calling riC{z, w) while using a global counter TV that records the number of nodes of the 
tree under construction. 

r/C(z, w): repeat 

N := 0; {counter for nodes} 

Can r!C{z,w) 

each time a node is built do 
N ■=N +1; 

if Bern((Ar + l)/(iV + 2)) continue; 

otherwise reject and restart from the first line; od 
until the generation finishes; 
return the object generated by T!C{z, w) 
(taking the distinguished leaf back into the set of U-atoms) 

Lemma 12. The algorithm TIC{z,w) is a Boltzmann sampler for the class K. of asym- 
metric bicolored binary trees. 

Proof. At each attempt, the call to TJC{z, w) would output a rooted binary tree 7 if there 
was no early interruption. Clearly, the probability that the generation of 7 finishes with- 
out interruption is YlfJi\i + ^)/ii + '2) = 2/{N{'y) + 2). Hence, each attempt is equivalent 
to doing 

7 <— ^!C{z, w); if Bern ^^^(^+2) I'^turn 7 else reject; 
Thus, the algorithm T]C(z,w) is equivalent to the algorithm given in the discussion pre- 
ceding Lemma [T^ hence r/C(0, w) is a Boltzmann sampler for the family IC. □ 

4.1.7. Boltzmann sampler for irreducible dissections. As stated in Proposition 1111 the 
closure-mapping realises a combinatorial isomorphism between asymmetric bicolored bi- 
nary trees (class K.) and asymmetric bicolored irreducible dissections (class T). Hence, 
the algorithm 

T2{z,w): T<— r/C(z,w); 

return closure(T) 

is a Boltzmann sampler for 2. In turn this easily yields a Boltzmann sampler for the 
corresponding rooted class JT. Precisely, starting from an asymmetric bicolored irreducible 
dissection, each of the 3 outer black vertices, which are distinguishable, might be chosen 
as the root-vertex in order to obtain a rooted irreducible dissection. Moreover the sets 
of L-atoms and U-atoms are slightly different for the classes X and JT"; indeed, a rooted 
dissection has one more L-atom (the black vertex following the root-vertex in cw order 
around the outer face) and one more U-atom (all faces are U-atoms in J', whereas only 
the inner faces are U-atoms in xj^. This yields the identity 

(10) J = 3*Zl*Zu*I, 

which directly yields (by the sampling rules of Figure [3]) a Boltzmann sampler TJ^{z,w) 
for from the Boltzmann sampler TX[z, w). 



We have chosen to specify the sets of L-atoms and U-atoms in this way in order to state the isomor- 
phisms /C ~ X and J7a — AI3. 
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Finally, we obtain a Boltzmann sampler for rooted admissible dissections by a simple 
rejection procedure 

TJi^{z, w): repeat 5 ^ TJ{z, w) until S G Ja', 
return S 



4.1.8. Boltzmann sampler for rooted 3-connected maps. The Boltzmann sampler for rooted 
irreducible dissections and the primal-map construction yield the following sampler for 
rooted 3-connected maps: 

rMl{z,w): S^rj,{z,w); 

return Primal((5) 

where Primal((5) is the rooted 3-connected map associated to S (see Section [4. 1.3|) . 

4.1.9. Boltzmann sampler for edge-rooted 3-connected planar graphs. To conclude, the 
Boltzmann sampler rA^3(z, w) yields a Boltzmann sampler TQ3(z,w) for edge-rooted 3- 
connected planar graphs, according to the isomorphism (Whitney's theorem) AI3 ~ 2*C/3, 

TQ3{z,'w): return rA^3(z, w) (forgetting the embedding) 

4.2. Boltzmann sampler for 2-connected planar graphs. The next step is to realise 
a Boltzmann sampler for 2-connected planar graphs from the Boltzmann sampler for edge- 
rooted 3-connected planar graphs obtained in Section 14.11 Precisely, we first describe 
a Boltzmann sampler for the class Q2 of edge-rooted 2-connected planar graphs, and 
subsequently obtain, by using rejection techniques, a Boltzmann sampler for the class G2 
of derived 2-connected planar graphs (having a Boltzmann sampler for Q2 allows us to go 
subsequently to connected planar graphs). 

To generate edge-rooted 2-connected planar graphs, we use a well-known decompo- 
sition, due to Trakhtenbrot [3T], which ensures that an edge-rooted 2-connected planar 
graph can be assembled from edge-rooted 3-connected planar components. This decom- 
position deals with so-called networks (following the terminology of Walsh [34]), where 
a network is defined as a connected graph N with two distinguished vertices and 00 
called poles, such that the graph N* obtained by adding an edge between and 00 is 
a 2-connected planar graph. Accordingly, we refer to Trakhtenbrot 's decomposition as 
the network decomposition. Notice that networks are closely related to edge-rooted 2- 
connected planar graphs, though not completely equivalent (see Equation below for 
the precise relation). 

We rely on |34| for the description of the network decomposition. A series-network 
or s-network is a network made of at least 2 networks connected in chain at their poles, 
the 00-pole of a network coinciding with the 0-pole of the following network in the chain. 
A parallel network or p-network is a network made of at least 2 networks connected in 
parallel, so that their respective 00-poles and 0-poles coincide. A pseudo-brick is a network 
TV whose poles are not adjacent and such that N* is a 3-connected planar graph with at 
least 4 vertices. A polyhedral network or /i-network is a network obtained by taking a 
pseudo-brick and substituting each edge e of the pseudo-brick by a network (polyhedral 
networks establish a link between 2-connectcd and 3-connccted planar graphs). 

Proposition 13 (Trakhtenbrot). Networks with at least 2 edges are partitioned into s- 
networks, p-networks and h-networks. 
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Let us explain how to obtain a recursive decomposition involving the different families 
of networks. (We simply adapt the decomposition formalised by Walsh [34] so as to 
have only positive signs.) Let V, S, V, and H be respectively the classes of networks, 
s-networks, p-networks, and /i-networks, where the L-atoms are the vertices except the 
two poles, and the U-atoms are the edges. In particular, Zjj stands here for the class 
containing the link-graph as only object, i.e., the graph with one edge connecting the two 
poles. Proposition [T3l ensures that 

V = Zu + S + V + H. 

An s-network can be uniquely decomposed into a non-s-network (the head of the chain) 
followed by a network (the trail of the chain) . which yields 



S = [Zu +v + n)*ZL*v. 

A p-network has a unique maximal parallel decomposition into a collection of at least 
two components that are not p-networks. Observe that we consider here graphs without 
multiple edges, so that at most one of these components is an edge. Whether there is one 
or no such edge-component yields 



V^Zu* Set>i(5 + n) + Set>2(5 4- H). 

By definition, the class of /i-networks corresponds to a U-substitution of networks in 
pseudo-bricks; and pseudo-bricks are exactly edge-rooted 3-connected planar graphs. As 
a consequence (recall that Ga stands for the family of 3-connected planar graphs), 



To sum up, we have the following grammar corresponding to the decomposition of 
networks into edge-rooted 3-connected planar graphs: 



(N) 



V = Zu + s + v + n, 

s = {Zu + v + n)*ZLi.v, 

V = Zu* Set>i(5 + n) + Set>2(5 

n = GlouV. 



W), 



Using the sampling rules (Figure [3]), the decomposition grammar (N) is directly trans- 
lated into a Boltzmann sampler T'D{z,y) for networks, as given in Figure [HI A network 
generated by TV{z, y) is made of a series-parallel backbone (3 (resulting from the branch- 
ing structures of the calls to rS{z, y) and rP(z, y)) and a collection of rooted 3-connected 
planar graphs that are attached at edges of f3; clearly all these 3-connected components are 
obtained from independent calls to the Boltzmann sampler TGaiz^w), with w = D{z,y). 

The only terminal nodes of the decomposition grammar are the classes Zl, Zij (which 
are explicit), and the class C/3. Thus, the sampler ri?(z,?/) and the auxiliary samplers 
r5(0,y), TV{z,y), and TT-L{z,y) are recursively specified in terms of TG'i{z,w), where w 
and z are linked by w = D[z, y). 

Observe that each edge-rooted 2-connected planar graph different from the link-graph 
gives rise to two networks, obtained respectively by keeping or deleting the root-edge. 
This yields the identity 



(11) 



{l+Zu)irG2 = (1+2?). 
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TV{z, y): Call TZu{z, y) or TS{z, y) or TV{z, y) or Vn{z, y) 
with respective probabilities -p, -^j ^, ^; 
return the network generated 

TZu{z^y): return the link-graph 

TS{z,y): ^ T{Zu +V + 'H){z,y); 
72 ^ ^(z,^); 
7 ^ 7i in series with 72; 
return 7 

r7'(z,y): Call VVi{z,y) or 

with resp. probabilities ^™p>^('^+^) ^ exp>2^s+£/). 
return the network generated 

TVi{z,y): fc ^ Pois>i(5 + i?); 

7i^r(5 + 7^)(z,w),...,7fe^r(5 + 7i)(z,i(;); {ind. calls} 

7 ^ (71, . . . , 7fc) in parallel; 

add to 7 an edge connecting the 2 poles; 

return 7 

fc^Pois>2(5 + ff); 

7i^r(5 + H)(2,w),...,7fe^r(5 + H)(z,w); {ind. calls} 
7 ^ (71, ... , 7fc) in parallel; 
return 7 

r7i:(z,y): 7 ^ TQziz^w), with tii = D{z,y); 
for each edge e of 7 do 
7e <- ^(z,?/); 
substitute e by 7e; 

{the poles of 7e are identified with the ends of e 
in a canonical way} 

od; 

return 7 

r{S + n)iz,y): Call TS{z,y) orrH(z,2/) 

with resp. probabilities ^ ^ ■ 



return the network generated 

r{Zu+V+n){z,y): Call rZu{z,y) OYrV{z,y) oi-Tn{z,y) 
with resp. probabilities 



H 



y+P+H' y+P+H^ y+P+H ' 

return the network generated 



Figure 9. Boltzmann samplers for networks. All generating functions 
are assumed to be evaluated at {z,y), i.e., D := D{z,y), S :~ S{z,y), 
P -Piz^y), and H:=H{z,y). 
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From that point, a Boltzmann sampler is easily obtained for the family Q2 of edge-rooted 
2-connected planar graphs. Define a procedure AddRootEdge that adds an edge con- 
necting the two poles and 00 of a network if they are not already adjacent, and roots 
the obtained graph at the edge (0,oo) directed from to 00. The following sampler for 
Q2 is the counterpart of Equation pT|) . 

r(l + 'D)(z, y): if Bern ^ i^ij{^z y) ) I'cturn the link-graph else return VDi^z, y); 
^G2{z, y): 7 ^ r(l + 'D){z, y); AddRootEdge(7); return 7 

Lemma 14. The algorithm TQ2{z,y) is a Boltzmann sampler for the class Q2 of edge- 
rooted 2-connected planar graphs. 

Proof. Firstly, observe that TQ2{z, y) outputs the link-graph either if the initial Bernoulli 
choice X is 0, or if X = 1 and the sampler VD{z,y) picks up the link-graph. Hence 
the link-graph is returned with probability (1 + y)/(l + D{z,y)), i.e., with probability 
1/G^(z,y). 

Apart from the link-graph, each graph j G G2 appears twice in the class £ := 1 -\- 
V: once in f |^| (keeping the root-edge) and once in (deleting the root- 

edge). Therefore, 7 has probability E{z,y)~^z^''^ /\^\\{y^^''^^'^^ -f j/"'''") of being drawn by 
TG2{z,y), where E{z,y) = 1 -f D{z,y) is the series of £. This probability simplifies to 
z^''^ /\j\\y^^^^^ /G2{z,y). Hence, rQ2{z,y) is a Boltzmann sampler for the class 02- □ 

The last step is to obtain a Boltzmann sampler for derived 2-connected planar graphs 
(i.e., with a distinguished vertex that is not labelled and does not count for the L-size) 
from the Boltzmann sampler for edge-rooted 2-connected planar graphs (as we will see 
in Section [4.31 derived 2-connected planar graphs constitute the blocks to construct con- 
nected planar graphs). 

We proceed in two steps. Firstly, we obtain a Boltzmann sampler for the U-derived 
class G2 (i-C, with a distinguished undirected edge that does not count in the U-size). 

Note that :— 2-kQ2 satisfies J- = ZL-kQ2. Hence, TQ2{z,y) directly yields a Boltzmann 
sampler TJ^{z,y) (sec the sampling rules in Figure [3]). Since T ~ 2-kQ2, a Boltzmann 
sampler for Q2 is obtained by calling TJ-(z, y) and then forgetting the direction of the 
root. 

Secondly, once we have a Boltzmann sampler TQ2{z^y) for the U-derived class Q2^ we 
just have to apply the procedure Uderived— >Lderived (described in Section [3.4. Sp to 
the class Q2 in order to obtain a Boltzmann sampler TQ2{z,y) for the L-derived class 
Q2. The procedure Uderived^Lderived can be successfully applied, because the ratio 
vertices/edges is bounded. Indeed, each connected graph 7 satisfies I7I < II7II -t~ 1, which 
easily yields a^/c/ = 2 for the class Q2 (attained by the link-graph). 

4.3. Boltzmann sampler for connected planar graphs. Another well known graph 
decomposition, called the block-decomposition, ensures that a connected graph can be 
decomposed into 2-connected components. We take advantage of this decomposition in 
order to specify a Boltzmann sampler for derived connected planar graphs from the Boltz- 
mann sampler for derived 2-connected planar graphs obtained in the last section. Then, 
a further rejection step yields a Boltzmann sampler for connected planar graphs. 

The block- decomposition (see [HJ p. 10] for a detailed description) ensures that each 
derived connected planar graph can be uniquely constructed in the following way: take a 
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set of derived 2-connected planar graphs and attach them together, by merging their 
marked vertices into a unique marked vertex. Then, for each unmarked vertex v of 
each 2-connected component, take a derived connected planar graph and merge the 
marked vertex of 7t, with v (this operation corresponds to an L-substitution). The block- 
decomposition gives rise to the following identity relating the classes Qi and Q2: 

(12) Sl' = SET(S2'oL(2L*gi')). 

This is directly translated into the following Boltzmann sampler for Q{ using the sampling 
rules of Figure [3l (Notice that the 2-connected blocks of a connected graph are built 
independently, each block resulting from a call to the Boltzmann sampler TQ-iiz, y), where 
z = xG{{x,y).) 

Tg{{x,y): k^Vois{Gd{z,y)); [With z = xG{{x,y)] 

7 ^ {TQ-iiz, ?/),..., TQ-iiz, y)); {k independent calls} 
merge the k components of 7 at their marked vertices; 
for each unmarked vertex v of 7 do 
Iv ^ TQ{{x,y)] 

merge the marked vertex of 7„ with v 

od; 

return 7. 

Then, a Boltzmann sampler for connected planar graphs is simply obtained from 
VQi{x, y) by using a rejection step so as to adjust the probability distribution: 

^Qi{x,y). repeat 7 <— rC/i(a::, y) 

take the marked vertex v back to the set of L-atoms; 
(if we consider the labels, v receives label I7I -I- 1) 
{this makes I7I increase by 1, and 7 G Qi} 

until Bern f ■; — r ) ; 

VI7I/ 

return 7 

Lemma 15. The sampler TQi{x^y) is a Boltzmann sampler for connected planar graphs. 

Proof. The proof is similar to the proof of Lemma [5] Due to the general property that 
Cn,m identifies to C'^_i ^, the sampler delimited inside the repeat/until loop draws each 

object J G Qi with probability G{(x, y)~^ (^^jli^i j/^^'^^L i-C-, with probability proportional 

to |7|^!yjT2/"''"- Hence, according to LemmalHl the sampler VQi{x,w) draws each object 

7 e CJi with probability proportional to jr^y"'*''!, i.e., is a Boltzmann sampler for Qi. □ 

4.4. Boltzmann sampler for planar graphs. A planar graph is classically decomposed 
into the set of its connected components, yielding 

(13) g = SET(gi), 

which translates to the following Boltzmann sampler for the class Q of planar graphs (the 
Set construction gives rise to a Poisson law, see Figure [S]): 

Tg{x,y): k^Vo\s{Gi{x,y))- 

return {TQi{x, y), . . . , rQi{x, y)) {k independent calls} 
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Proposition 16. The procedure TQ{x,y) is a Boltzmann sampler for planar graphs. 

5. Deriving an efficient sampler 

We have completely described in Section [5] a mixed Boltzmann sampler rC/(a;,y) for 
planar graphs. This sampler yields an exact-size uniform sampler and an approximate-size 
uniform sampler for planar graphs: to sample at size n, call the sampler TQ{x^ 1) until 
the graph generated has size n; to sample in a range of sizes [n(l — e),n(l -I- e)], call the 
sampler VQ{x, 1) until the graph generated has size in the range. These targetted samplers 
can be shown to have expected polynomial complexity, of order n^/'^ for approximate-size 
sampling and rJ/"^ for exact-size sampling (we omit the proof since wc will describe more 
efficient samplers in this section). 

However, more is needed to achieve the complexity stated in Theorem [1] i.e., 0(n/e) 
for approximate-size sampling and 0{n'^) for exact-size sampling. The main problem of 
the sampler TQix^ 1) is that the typical size of a graph generated is small, so that the 
number of attempts to reach a large target size is prohibitive. 

In order to correct this effect, we design in this section a Boltzmann sampler for "bi- 
derived" planar graphs, which are equivalent to bi-pointed planar graphs, i.e., with 2 
distinguished verticefl The intuition is that a Boltzmann sampler for bi-pointed planar 
graphs gives more weight to large graphs, because a graph of size n gives rise to n{n — 1) 
bi-pointcd graphs. Hence, the probability of reaching a large size is better (upon choosing 
suitably the value of the Boltzmann parameter). The fact that the graphs have to be 
pointed 2 times is due to the specific asymptotic behaviour of the coefficients counting 
planar graphs, which has been recently analysed by Gimenez and Noy [20j . 

5.1. Targetted samplers for classes with square-root singularities. As we describe 
here, a mixed class C with a certain type of singularities (square-root type) gives rise to 
efficient approximate-size and exact-size samplers, provided C has a Boltzmann sampler 
such that the expected cost of generation is of the same order as the expected size of the 
object generated. 

Definition 17. Given a mixed class C, we define a singular point of C as a pair xq > 0, 
2/0 > such that the function x i-^ C(.t,?/o) has a dominant singularity at xq (the radius 
of convergence is xq). 

Definition 18. For a e M\Z>o, a mixed class C is called a-singular if, for each singular 
point (xo,yo) of C, the function a; i-^ C{x,yo) has a unique dominant singularity at xq 
(i.e., Xq is the unique singularity on the circle \z\ = xq) and admits a singular expansion 
of the form 

C{x, yo) = P{x) + Ca ■ {xq -x)" +o ((xo - x)") , 

where Cq, is a constant, P{x) is rational with no poles in the disk \z\ < xq, and where the 
expansion holds in a so-called A-neighbourhood of xq, see [14[ [T5] . In the special case 
a ~ 1/2, the class is said to have square-root singularities. 

^In an earlier version of the article and in the conference version |16| , we derived 3 times — as prescribed 
by 111! — in order to get a singularity type (1 — x/p)~^/^ (efficient targetted samplers are obtained when 
taking x = p(f — f/(2n))). We have recently discovered that deriving 2 times (which yields a square-root 
singularity type (f — x/p)^/^) and taking again x = p{l — l/{2n)) yields the same complexities for the 
targetted samplers, with the advantage that the description and analysis is significantly simpler (in the 
original article |lf) . they prescribe to take x = p and to use some early abort techniques for square-root 
singularity type, but it seems difficult to analyse the gain due to early abortion here, since the Boltzmann 
sampler for planar graphs makes use of rejection techniques). 



UNIFORM RANDOM SAMPLING OF PLANAR GRAPHS IN LINEAR TIME 



29 



Lemma 19. Let C be a mixed class with square-root singularities, and endowed with a 
Boltzmann sampler rC(x, y). Let (xo, yo) be a singular point of C. For any n > Q, define 



Call 7r„ resp.) the probability that an object 7 generated by TC{xn,yo) satisfies 

\l\ = (I7I S := — e), n(l + e)], resp.); and call an the expected size of the output 
ofrC{xn,yo)- 

Then l/7r„ is 0{n^^^), l/7r„^g is 0{n^^'^ /e), and On is 

Proof. The so-called transfer theorems of singularity analysis [13] ensure that the coeffi- 
cient a„ := [a;"]C(a;, j/o) satisfies, as n ^ 00, a„ ~„^oo c where c is a positive 
constant. This easily yields the asymptotic bounds for l/iTn and l/7r„_e, using the expres- 
sions 7r„ = a„x„"/C(a;„,?/o) and 7r„^e = J2kei„,, akxj' /C{x„,yo). 

It is also an easy exercise to find the asymptotics of (t„, using the formula (given in 

■ dxC{xn,yo)/C{xn,yo)- □ 

Lemma [T51 suggests the following simple heuristic to obtain efficient targetted samplers. 
For approximate-size sampling (exact-size sampling, resp.), repeat calling TC{xn, 1) until 
the size of the object is in (is exactly n, resp.). (The parameter y is useful if a target 
U-sizc m is also given, as we will see for planar graphs in Section 16.21 ) The complexity 
of sampling will be good for a class C that has square-root singularities and that has 
an efficient Boltzmann sampler. Indeed, for approximate-size sampling, the number of 
attempts to reach the target-domain /„ (i.e., t:~\) is of order n^/^, and for exact-size 
sampling, the number of attempts to reach the size n (i.e., tt"^) is of order n"^/^. If 
C is endowed with a Boltzmann sampler TC{x,y) such that the expected complexity of 
sampling at {xn,yo) is of order y/n (same order as the expected size ct„), then the expected 
complexity is typically 0{n/e) for approximate-size sampling and 0{n^) for exact-size 
sampling, as we will see for planar graphs. 

Let us mention that the original article [11] uses a different heuristic. The targetted 
samplers also repeat calling the Boltzmann sampler until the size of the object is in the 
target domain, but the parameter x is chosen to be exactly at the singularity p. The 
second difference is that, at each attempt, the generation is interrupted if the size of the 
object goes beyond the target domain. We prefer to use the simple heuristic discussed 
above, which does not require early interruption techniques. In this way the samplers are 
easier to describe and to analyse. 

In order to apply these techniques to planar graphs, we have to derive two times the 
class of planar graphs, as indicated by the following two lemmas. 

Lemma 20 ( [1A\ ) . If a class C is a-singular, then the class C is (a — \)-singular (by the 
effect of derivation) . 

Lemma 21 ([20j). The class Q of planar graphs is b/2-singular, hence the class Q" of 
bi-derived planar graphs has square-root singularities. 

5.2. Derivation rules for Boltzmann samplers. As suggested by Lemma [12] and 
Lemma I21[ we will get good targetted samplers for planar graphs if we can describe an 
efficient Boltzmann sampler for the class G" of bi-derived planar graphs (a graph in Q" 
has two unlabcUed vertices that are marked specifically, say the first one is marked * and 
the second one is marked Our Boltzmann sampler TQ"{x, y) — to be presented in this 
section — makes use of the decomposition of planar graphs into 3-connected components 
which we have already successfully used to obtain a Boltzmann sampler for planar graphs 
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in SectionlH This decomposition can be formally translated into a decomposition grammar 
(with additional unpointing/pointing operations). To obtain a Boltzmann sampler for bi- 
derivcd planar graphs instead of planar graphs, the idea is simply to derive this grammar 
2 times. 

As we explain here and as is well known in general, a decomposition grammar can 
be derived automatically. (In our framework, a decomposition grammar involves the 5 
constructions {+,*, SET>rf, o^, oy}.) 

Proposition 22 (derivation rules). The basic finite classes satisfy 

(1)' = 0, (Zi)' = l, {Zuy = 0. 

The 5 constructions satisfy the following derivation rules: 

(A + BY = A' + B', 

{A*B)' = A'i^B + Ai^B', 

(14) { (SET>d(S))' = S'*SET>rf_i(6) for d> 0, (with Set>_i = Set) 

{AolB)' = 6'*(^'oiS), 

{AouB)' ^ A' oijB + B' i.{AouB). 

Proof. The derivation formulas for basic classes are trivial. The proof of the derivation 
rules for {+,*, o^} are given in [3]. Notice that the rule for SET>ci follows from the rule 
for ol. (Indeed, SET>£i(S) = Aoj^ B, where A = SET>d(Zi), which clearly satisfies 
A' = SET>d_i (Zi).) Finally, the proof of the rule for uses similar arguments as the 
proof of the rule for oj^. In an object of {Aojj B)' , the distinguished atom is either on the 
core-structure (in A), or is in a certain component (in B) that is substituted at a certain 
U-atom of the core-structure. The first case yields the term A' ojj B, and the second case 
yields the term B' *{Aou B). □ 

According to Proposition [22l it is completely automatic to find a decomposition gram- 
mar for a derived class C if we are given a decomposition grammar for C. 

5.3. Boltzmann sampler for bi-derived planar graphs. We present in this section 
our Boltzmann sampler TQ"{x,y) for bi-derived planar graphs, with a quite similar ap- 
proach to the one adopted in Section [4l and again a bottom-to-top presentation. At 
first the closure-mapping allows us to obtain Boltzmann samplers for 3-connected planar 
graphs marked in various ways. Then we go from 3-connected to bi-derived planar graphs 
via networks, bi-derived 2-connected, and bi-derived connected planar graphs. 

The complete scheme is illustrated in Figure [TUl which is the counterpart of Figure [H 

5.3.1. Boltzmann samplers for derived binary trees. We have already obtained in Sec- 
tion |4T3] a Boltzmann sampler for the class /C of unrooted asymmetric binary trees. Our 
purpose here is to derive a Boltzmann sampler for the derived class K.'. Recall that we 
have also described in Section [4. 1.51 a Boltzmann sampler for the U-derived class IC, which 
satisfies the completely recursive decomposition grammar ([9]) (see also Figure H]). Hence, 
we have to apply the procedure Uderived-^Lderived described in Section r3. 4. 31 to the 
class /C in order to obtain a Boltzmann sampler TIC'{z,w) from TIC{z,w). For this we 
have to check that aj^/u is finite for the class K,. It is easily proved that a bicolored binary 
tree with m leaves has m — 2 nodes, and that at most [2(m — 3)/3j of the nodes are black. 
In addition, there exist trees with 3i -I- 3 leaves and 2i black nodes (those with all leaves 
incident to black nodes). Hence, for the class /C, the parameter aL/jj is equal to 2/3. 
Therefore the procedure Uderived^Lderived can be apphed to the class /C. 
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Figure 10. The complete scheme to obtain a Boltzmann sampler for 
bi-dcrived planar graphs. 



5.3.2. Boltzmann samplers for derived rooted dissections and 3-connected maps. Our next 
step is to obtain Boltzmann samplers for derived irreducible dissections, in order to go 
subsequently to 3-conncctcd maps. As expected we take advantage of the closure-mapping. 
Recall that the closure-mapping realises the isomorphism /C ~ X between the class /C of 
asymmetric binary trees and the class I of asymmetric irreducible dissections. There is no 
problem in deriving an isomorphism, so the closure-mapping also realises the isomorphism 
/C' ~ X'. Accordingly we have the following Boltzmann sampler for the class I': 

T2'{z,w): T^TIC'{z,w); 

6 ^ closure(T); 
return S 

where the discarded L-atom is the same in r and in S. 

Then, we easily obtain a Boltzmann sampler for the corresponding rooted class J^'. 
Indeed, the equation J = 3-kZL* Zu *I that relates X and J yields J' = 3* Zu *I -I- 3 ★ 
ZL*ZifkT' . Hence, using the sampling rules of Figure [3l we obtain a Boltzmann sampler 
TJ''(z,w) from the Boltzmann samplers TX[z,w) and TX'(z,w). 
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From that point, we obtain a Boltzmann sampler for the derived rooted dissections that 
are admissible. As J7a C X, we also have J!^ C J' , which yields the following Boltzmann 
sampler for J!^: 

Tj!^{z, w): repeat 6 ^ TJ'{z, w) 
until ^ e J'a; 
return 5 

Finally, using the isomorphism J7a — A^3 (primal map construction, Section r4.1.3p . 
which yields J!^ ~ A^s', we obtain a Boltzmann samplers for derived rooted 3-connected 
maps: 

return Primal((5) 

where the returned rooted 3-connected map inherits the distinguished L-atom of 6. 



5.3.3. Boltzmann samplers for derived rooted 3-connected planar graphs. As we have seen 
in Section I4.1.2[ Whitney's theorem states that any 3-connected planar graph has two 
embeddings on the sphere (which differ by a reflection). Clearly the same property holds 
for 3-connected planar graphs that have additional marks. (We have already used this 
observation in Section l4.1.2l for rooted graphs, M.^ ~ 2*C/3, in order to obtain a Boltzmann 
sampler for ^^3.) Hence A^s' ~ 2 ^^3', which yields the following Boltzmann sampler for 



Vg^'iz^w): return rAl3'(z,w); 

(forgetting the embedding) 



The next step (in Section is to go to derived networks. This asks for a derivation 
of the decomposition grammar for networks, which involves not only the classes Q^, Q^' , 
but also the U-dcrived class ^^3. Hence, we also need a Boltzmann sampler for Q^. 

To this aim wc just have to apply the procedure Lderived^Uderived to the class Q^- 
By the Eulcr relation, a 3-conncctcd planar graph with n vertices has at most 3n — 6 edges 
(equality holds for triangulations). Hence, the parameter ajz/^ is equal to 3 for the class 

Qz, so Lderived^Uderived can be successfully applied to Qz, yielding a Boltzmann 
sampler for from the Boltzmann sampler for Q^' . 

5.4. Boltzmann samplers for derived networks. Following the general scheme shown 
in Figure I10[ our aim is now to obtain a Boltzmann samplers for the class V of derived 
networks. Recall that the decomposition grammar for T) has allowed us to obtain a 
Boltzmann sampler for T) from a Boltzmann sampler for Q^. Using the derivation rules 
(Proposition injected in the grammar (N), we obtain the following decomposition 
grammar for P': 

-D' = s' + v' + n\ 

s' = {v + 'H')-kZ-kV + {c + v + n)i^{v + z-kV), 

V ^Ci^{S' + W) ★ Set(5 + H) + (5' + W) ★ Set>i(5 + U), 



(N') 
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The only terminal classes in this grammar arc Q3' and Gs- Hence, the sampling rules 
of Figure [3] yield a Boltzmann sampler for 2?' from the Boltzmann samplers for and 
Gs which wc have obtained in Section fS. 3. 31 The sampler n^'{z,y) looks similar (though 
with more cases) to the one for T'D{z, y) given in Figure [HI 

5.5. Boltzmann samplers for bi-derived 2-connected planar graphs. The aim of 

this section is to obtain Boltzmann samplers for the class G2 of bi-derivcd 2-connected 
planar graphs (after the Boltzmann sampler for G2 obtained in Section 14. 2[) , in order to 
go subsequently to bi-derivcd connected planar graphs. 

At first, the Boltzmann sampler for V' yields a Boltzmann sampler for the class G2 ■ 
Indeed the identity {l + V) = (1 + Zi/)-kG2 is derived as P' = (1 + Zij)-kG2', which yields 
the following sampler, 

r^'(z,y): j^TV'iz,y); 

AddRootEdge(7) ; 
return 7 

where AddRootEdge has been defined in Section B?^ The proof that this is a Boltzmann 
sampler for G2' is similar to the proof of Lemma [Ml 

Next we describe a Boltzmann sampler for the class G2'- As we have seen in Scction H?^ 

G2 and G2 are related by the identity 2-kG2 = *G2- Hence, if we define 2*^/2, we 
have T' = Zj} -kG2 + 2*Zi *C/2. Hence, the sampling rules of Figure [3] yield a Boltzmann 
sampler rjF'(z, y) for the class T' . Clearly, as JF' = 2 * G2 ^ a Boltzmann sampler for G2 
is obtained by calling VT'{z, y) and forgetting the direction of the root. 

Finally, the procedure Uderived— >Lderived yields (when applied to G2') from the 
Boltzmann sampler for G2' to a Boltzmann sampler for G2" ■ The procedure can be suc- 
cessfully applied, as the class G2 satisfies = 1 (attained by the link-graph). 

5.5.1. Boltzmann sampler for bi-derived connected planar graphs. The block-decomposition 
makes it easy to obtain a Boltzmann sampler for the class Gi" of bi-derived connected 
planar graphs (this decomposition has already allowed us to obtain a Boltzmann sampler 
for Gi in Section |473|) . Recall that the block-decomposition yields the identity 

Gi=Set{G^ol (Zl^Gi)), 

which is derived as 

= (G^ + Zl^ Gi') * G2 OL {Zl * Gi) * Gi. 

As we already have Boltzmann samplers for the classes G2" and Gi' , the sampling 
rules of Figure [3] yield a Boltzmann sampler rGi"{x,y) for the class Gi"- Observe that 
the 2-connectcd blocks of a graph generated by TGi"{x,y) are obtained as independent 
calls to rG2{z,y) and TG2" {z,y), where z and x are related by the change of variable 
z = xG{{x,y). 

5.5.2. Boltzmann samplers for bi-derived planar graphs. We can now achieve our goal, 
i.e., obtain a Boltzmann sampler for the class G" of bi-derivcd planar graphs. For this 
purpose, we simply derive twice the identity 

G = SET(gi), 

which yields successively the identities 

G' ^Gi^G, 
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and 

From the first identity and Vgix, y), Tg{{x, y), we get a Boltzmann sampler Tg'{x, y) for 
the class g' . Then, from the second identity and Vg{x, y), rg'{x, y), Tg{{x, y), rg{'{x, y), 
we get a Boltzmann sampler rg"{x, y) for the class g" . 

6. The targetted samplers for planar graphs 

The Boltzmann sampler Tg"{x,y) — when tuned as indicated in Lemma 1191 — yields 
efficient exact-size and approximate-size random samplers for planar graphs, with the 
complexities as stated in Theorem [T] and Theorem [2] Define the algorithm: 

SAMPLEPLANAR(a;, y): 7 ^ Vg"{x, y); 

give label I7I + 1 to the vertex marked * 
and label I7I + 2 to the marked vertex * 
(thus I7I increases by 2, and 7 G g)\ 
return 7 

6.1. Samplers according to the number of vertices. Let pc be the radius of con- 
vergence oi X ^ G{x,l). Define 

— (1 - h) ■ PG- 

For n > 1, the exact-size sampler is 

2l„: repeat 7 ^ SAMPLEPLANAR(a;„, 1) until I7I — n; return 7. 
For n > 1 and e > 0, the approximate-size sampler is 

2t„^e: repeat 7 ^ SamplePlanar(x„, 1) until I7I £ [?i(l - e),?i(l + e)]; return 7. 

6.2. Samplers according to the numbers of vertices and edges. For any y > 0, 

we denote by pciy) the radius of convergence of a; s- G{x,y). Let fi{y) be the function 
defined as 

Ky) -y'^iy)/PGiy)- 

dy 

As proved in (using the so-called quasi-power theorem), for a fixed y > 0, a large graph 
drawn by the Boltzmann sampler Tg"{x, y) has a ratio edges/ vertices concentrated around 
the value ^{y) as x approaches the radius of convergence of x 1-^ G{x,y). This yields a 
relation between the secondary parameter y and the ratio edges/vertices. If we want a 
ratio edges/ vertices close to a target value /i, we have to choose y so that p.{y) = p. It is 
shown in [20j that the function /i(j/) is strictly increasing on (0,-fcx)), with lini/i(j/) — 1 
as y — > and \im fi{y) = 3 as y ^ +00. As a consequence, fi{y) has an inverse function 
y(/z) defined on (1, 3). (In addition, fi 1-^ vifJ-) can be evaluated with good precision from 
the analytic equation it satisfies.) We define 

XnilJ-) (1 - ^) • PciyifJ-))- 
For n > 1 and fi G (1, 3), the exact-size sampler is 

Sln.^t! repeat 7 ^ SAMPLEPLANAR(2;„(/i), ?/(/i)) until (171=71. and | I7I | = Lmt^J ); return 7. 
For 71 > 1, ^ G (1, 3), and e > 0, the approximate-size sampler is 
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2ln,Ai,e: repeat 7 <— SamplePlanar(.t„(/^), ?;(^)) 

until (I7I e [n{l - e), n{l + e)] and ^ G [//(I - e),fi{l + e)]); 
return 7. 

The complexity of the samplers is analysed in Section [51 

7. Implementation and experimental results 

7.1. Implementation. We have completely implemented the random samplers for pla- 
nar graphs described in Section O First we evaluated with good precision — typically 20 
digits — the generating functions of the families of planar graphs that intervene in the de- 
composition (general, connected, 2-connectcd, 3-connected) , derived up to 2 times. The 
calculations have been carried out in Maple using the analytic expressions of Gimcnez 
and Noy for the generating functions [20]. We have performed the evaluations for values 
of the parameter x associated with a bunch of reference target sizes in logarithmic scale, 
n = {10^, 10^, 10"*, 10^, 10^}. From the evaluations of the generating functions, we have 
computed the vectors of real values that are associated to the random choices to be per- 
formed during the generation, e.g., a Poisson law vector with parameter Gi{x) (the EGF 
of connected planar graphs) is used for drawing the number of connected components of 
the graph. 

The second step has been the implementation of the random sampler in Java. To 
build the graph all along the generation process, it proves more convenient to manipulate 
a data structure specific to planar maps rather than planar graphs. The advantage is 
also that the graph to be generated will be equipped with an explicit (arbitrary) planar 
embedding. Thus if the graph generated is to be drawn in the plane, we do not need 
to call the rather involved algorithms for embedding a planar graph. Planar maps are 
suitably manipulated using the so-called half-edge structure, where each half-edge occupies 
a memory block containing a pointer to the opposite half-edge along the same edge and to 
the next half-edge in ccw order around the incident vertex. Using the half-edge structure, 
it proves very easy to implement in cost 0(1) all primitives used for building the graph — 
typically, merging two components at a common vertex or edge. Doing this, the actual 
complexity of implementation corresponds to the complexity of the random samplers as 
stated in Theorem [T] and Theorem [2| linear for approximate-size sampling and quadratic 
for exact-size sampling. In practice, generating a graph of size of order 10^ takes a few 
seconds on a standard computer. 

7.2. Experimentations. The good complexity of our random samplers allows us to ob- 
serve statistical properties of parameters on very large random planar graphs — in the 
range of sizes 10^ — where the asymptotic regime is already visible. We focus here on pa- 
rameters that are known or expected to be concentrated around a limit value. Note that 
the experimentations are on connected planar graphs instead of general planar graphs. 
(It is slightly easier to restrict the implementation to connected graphs, which are con- 
veniently manipulated using the half-edge data structure.) However, from the works of 
Gimenez and Noy [20j and previous work by MacDiarmid et al. [23j . a random planar 
graph consists of a huge connected component, plus other components whose total ex- 
pected size is 0(1). Thus, statistical properties like those stated in Conjecture E51 should 
be the same for random planar graphs as for random connected planar graphs. 

Number of edges. First we have checked that the random variable Xn that counts the 
number of edges in a random connected planar graph with n vertices is concentrated. 
Precisely, Gimenez and Noy have proved that F„ := Xn/n converges in law to a constant 
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Figure 11. Ratio edges/ vertices observed on a collection 71, . . . ,730 of 
80 random connected planar graphs of size at least lO'''; each graph 7^ 
yields a point at coordinates («,Rat(7i)), where Rat (7) is the ratio given 
by the number of edges divided by the number of vertices of 7. 



fi w 2.213, (they also show that the fluctuations are gaussian of magnitude 1/y/n). Fig- 
ured!] shows in ordinate the ratio edges/vertices for a collection of 80 random connected 
planar graphs of size at least 10^ drawn by our sampler. As we can see, the ratios are 
concentrated around the horizontal line y = fi, agreeing with the convergence result of 
Gimenez and Noy. 

Degrees of vertices. Another parameter of interest is the distribution of the degrees of 
vertices in a random planar graph. For a planar graph 7 with n vertices, we denote by 
]V'('=)(7) the number of vertices of 7 that have k neighbours. Accordingly, Z'^^\'^) := 
N^^\'^) /n is the proportion of vertices of degree k in 7. It is known from Gimenez 
and Noy that, for k = 1,2, the random variable Z'^''^ converges in law to an explicit 
constant. Figure [T^ shows in abscissa the parameter k and in ordinate the value of Z^*^^ 
for a collection of 80 random connected planar graphs of size at least lO** drawn by our 
sampler. Hence, the vertical line at abscissa k is occupied by 80 points whose ordinates 
correspond to the values taken by Z^''^ for each of the graphs. As we can see, for k 
small — typically k « logn — the values of Z^'^'f are concentrated around a constant. This 
leads us to the following conjecture. 

(k) 

Conjecture 23. For every k > 1, let Zn be the random variable denoting the proportion 
of vertices of degree k in a random planar graph with n vertices taken uniformly at random. 
Then zli^ converges in law to an explicit constant tt^'''^ as n —> 00; and tt'-'^-' = 1. 
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Figure 12. The distribution of vertex degrees observed on a 
collection 71 , . . . , 730 of 80 random connected planar graphs of 
size at least 10"*. Each graph 7 yields points at coordinates 
(l,Z(i)(7)),(2,Z(2)(7)),...,(d,Z('')(7)), where d is the maximal degree 
of 7 and, for 1 < fc < d, Z^''\'j) is the proportion of vertices of 7 that 
have degree k. 



Let us mention some progress on this conjecture. It has recently been proved in [10] 
that the expected values E(z4'^-') converge as 71 00 to constants tt'^'^'' that are computable 
and satisfy X]fe = 1- Hence, what remains to be shown regarding the conjecture is 
the concentration property. 



8. Analysis of the time complexity 

This whole section is dedicated to the proof of the complexities of the targettcd random 
samplers. We show that the expected complexities of the targetted samplers 2l„, 2l„_e, 
2t„^p, and 2l„_^_e, as described in Section[6l are respectively 0{n^), 0{n/e), Op(n^/^), and 
Ofj_{n/e) respectively (the dependency in fi in not analysed for the sake of simplicity). 

Recall that the targetted samplers call rQ"{x, y) (with suitable values of x and y) until 
the size parameters arc in the target domain. Accordingly, the complexity analysis is 
done in two steps. In the first step, we estimate the probability of hitting the target 
domain, which allows us to reduce the complexity analysis to the analysis of the expected 
complexity of the pure Boltzmann sampler TQ"{x, y). We use a specific notation to denote 
such an expected complexity: 
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Definition 24. Given a class C endowed with a Boltzmann sampler VC{x, y), we denote by 
AC(x, y) the expected combinatorial complexitjQ of a call to rC(x, y) (note that AC{x, y) 
depends not only on C, but also on a specific Boltzmann sampler for C). 

Typically the values {x,y) have to be close to a singular point of Q in order to draw 
graphs of large size. Hence, in the second step, our aim is to bound KQ" [x^ y) when (x, y) 
converges to a given singular point (a;o,?/o) of Q. To analyse KQ"{x,y), our approach 
is again from bottom to top, as the description of the sampler in Section [5] (see also 
the general scheme summarized in Figure llOp . At each step we give asymptotic bounds 
for the expected complexities of the Boltzmann samplers when the parameters approach 
a singular point. This study requires the knowledge of the singular behaviours of all 
series involved in the decomposition of bi-derived planar graphs, which are recalled in 
Section 

8.1. Complexity of rejection: the key lemma. The following simple lemma will be 
extensively used, firstly to reduce the complexity analysis of the targetted samplers to the 
one of pure Boltzmann samplers, secondly to estimate the effect of the rejection steps on 
the expected complexities of the Boltzmann samplers. 

Lemma 25 (rejection complexity). Let 2t he a random sampler on a combinatorial class 
C according to a probability distribution P, and let p : C be a function on C , called 

the rejection function. Consider the rejection algorithm 
Slrcj- repeat 7 <— 21 until Bern(p(7)) return 7. 

Then the expected complexity E(2lroj) o/2lioj and the expected complexity E(2l) o/2l are 
related by 

(15) E(2l,ej) = — E(a), 

Pace 

where p^cc S7GC ^(7)^(7) '■^ probability of success o/2lrcj at each attempt. 
Proof. The quantity E(2trcj) satisfies the recursive equation 

E(2l,ej) = E(2l) + (1 - Pacc)E(2l,.cj). 

Indeed, a first attempt, with expected complexity E(2t), is always needed; and in case of 
rejection, occurring with probability (1 — Pacc), the sampler restarts in the same way as 
when it is launched. □ 

As a corollary we obtain the following useful formulas to estimate the effect of rejection 
in Boltzmann samplers when going from L-derived (vertex-pointed) to U-derived (edge- 
pointed) graphs and vice-versa. 

Corollary 26 (Complexity of changing the root). Let A be a mixed combinatorial class 
such that the constants ctu/i^ := max^g_4-^|^ and aj^ju := max-y£_4y|^ are finite. Define 
c := au/L ■ OiL/u- 

• Assume A' is equipped with a Boltzmann sampler, and let TA{x, y) be the Boltz- 
mann sampler for A obtained by applying Lderived^Uderived — as defined in 
Section \3jl\ —to A. Then 

A-4(a;,y) < c • AA'{x,y). 



'See the discussion on the complexity model after the statement of Theorem [2] in the introduction. 
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• Assume A is equipped with a Boltzmann sampler, and let VA!{x,ii) he the Boltz- 
mann sampler for A' obtained by applying U DERIVED—!- L DERIVED — as defined 
in Section \S.4-^ — to A. Then 

AA'{x,y) < c- AA{x,y). 

Proof. Let us give the proof for Lderived^Uderived (the other case is proved in a 
similar way). By definition of Lderived^Uderived the probabihty of the BernouUi 
choice at each attempt mTA{x,y) is at least min^g_4Jj^, i.e., at least l/{au/L''^L/u)- 
Hence the probability Pacc of success at each attempt is at least 1/c. Therefore, by 
CoroUaryHi AA{x, y) ^ AA!{x, y)/pacc < c • AA!{x, y). □ 

8.2. Reduction to analysing the expected complexity of Boltzmann samplers. 

We prove here that analysing the expected complexities of the targetted samplers reduces 
to analysing the expected complexity AG"{x,y) when {x,y) approaches a singular point. 
(Recall that a singular point (xo,yo) for a class C is such that the function x t—f C{x,yo) 
has a dominant singularity at Xo-) 

Claim 27. Assume that for every singular point {xo,yo) ofQ, the expected complexity of 
the Boltzmann sampler for Q" satisfied 

(16) Ae"(x,?/o) = 0((a;o-x)-i/2) asx->xo. 

Then the expected complexities of the targetted samplers 2l„, 2l„_e, 2ln,/j, and ^n,fj.,e — as 
defined in Section\^ — are respectively 0{n'^), 0{n/e), O^in^/"^), and 0^{n/e). 

In other words, proving p6p is enough to prove the complexities of the random samplers 
for planar graphs, as stated in Theorem [7] and Theorem O 

Proof. Assume that (|T6l) holds. Let 7r„^e (TTn, resp.) be the probability that the output of 
SAMPLEPLANAR(a;„, 1) — with Xn = (1 — l/2n)-pG — has size in J„_e := [n(l — e), n(l + e)] 
(has size n, resp.). According to Lemma the expected complexities of the exact-size 
and approximate-size samplers with respect to vertices — as described in Section 16.11 — 
satisfy 

Equation (|16p ensures that, when n oo, AQ"{xn, 1) is 0{n^^^). In addition, according 
to Lemma im G" is 1/2-singular (square- root singularities). Hence, by Lemma [TOl l/7r„ 
is 0(n3/2) ^j^^^^ ig 0(ni/Ve). Thus, E(an) is 0(j?) and E(2ln,e) is (^(n/e). 

The proof for the samplers with respect to vertices and edges is a bit more technical. 
Consider a planar graph 7 drawn by the sampler SAMPLEPLANAR(a;„(/i), ?/(/i)). In view 
of the proof for the exact-size sampler, define 

^uAM :=]P(II7||=LmH,|7| = "), Tf^in :=P(||7lhLA^nJ \ \l\^n), 7r„ P(|7|=n). 
In view of the proof for the approximate-size sampler, define 

^„AM,. P(l7l e [^(1 - e),"(l + e)], II7II/I7I e [/i(l - e),/i(l + e)]), 

^Min,. ■■= miWIH e [m(i - e), m(i + e)] I l7l e [«(i - e), ^(1 + e)]), 

and 

7r„., :=P(|7| e [n(l-e),n(l + e)]). 



'^In this article all convergence statements are meant "from below", i.e., a; — » means that x ap- 
proaches XQ while staying smaller than xq. 
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Notice that TTnAfi ~ ^^fiin ■ and nnA^.e = TTpin^e ' 7r„^e. Morcovcr, Lemma [55] ensures that 

It has been shown by Gimenez and Noy [20] (based on the quasi-power theorem) that, 
for a fixed fi £ (1, 3), l/7f^|„ is Op(n^/^) as ^ oo (the dependency in /i is not discussed 
here for the sake of simphcity). Moreover, Lemma [TOl ensures that l/7r„ is 0^(n^/^) as 
n oo. Hence, l/7f„^^ is Of^{n^). FinaUy Equation ([IB)) ensures that At/"(a;„(/Lt), j/(/x)) 
is Op(ni/2), therefore E(2l„,^) is ©^(n^/^). 

For the approximate-size samplers, the results of Gimenez and Noy (central limit theo- 
rems) ensure that, when /i £ (1,3) and e > are fixed and n oo, 7fp|„ ^ converges to 1. 
In addition. Lemma [T9l ensures that l/nn^^ is 0^(?i^/^/e). Hence, l/7f„A/^,c is Ofj,{n^^^ /e). 
Equation implies that ACJ"(x„(/i), ?/(^)) is ©^(n^/^), hence E(2l„,p,e) is 0^,{n/e). □ 

From now on, our aim is to prove that, for any singular point (xq, ?/o) of Qi AC/"(x, yo) 
is 0{{xq — x)^^^"^) as X ^ .To- 

8.3. Expected sizes of Boltzmann samplers. Similarly as for the expected complex- 
ities, it proves convenient to use specific notations for the expected sizes associated to 
Boltzmann samplers, and to state some of their basic properties. 

Definition 28 (expected sizes). Let C be a mixed combinatorial class, and let {x,y) be 
admissible for C (i.e., C{x,y) converges). Define respectively the expected L-size and the 
expected U-sizc at (x, y) as the quantities 



'^L.l|7|l_^5-C'(a;,y) 



a;- 



C(x,y)^'"|7r C{x,y) ' 

We will need the following two simple lemmas at some points of the analysis. 

Lemma 29 (monotonicity of expected sizes). Let C be a mixed class. 

• For each fixed yo > 0, the expected L-size x |C|(2, is increasing with x. 

• For each fixed xq > 0, the expected U-size y > |C|(3;g j,) is increasing with y. 

Proof. As noticed in jllj (in the labelled framework), the derivative of the function f{x) := 
\C\{x,yo) is equal to 1/x multiplied by the variance of the L-size of an object under the 
Boltzmann distribution at {x^yo). Hence f'{x) > for x > 0, so f{x) is increasing with 
X. Similarly the derivative of g{y) :~ \\C\\(xo,y) is equal to 1/y multiplied by the variance 
of the U-size of an object under the Boltzmann distribution at (xo,y), hence g{y) is 
increasing with y for y > 0. □ 

Lemma 30 (divergence of expected sizes at singular points). Let C be an a-singular class 
and let (xo, j/o) be a singular point of C. Then, as x ^ xq.' 

u if a > 1, the expected size x ^ \C\{x,yo) converges to a positive constant, 

• i/0 < a < 1, the expected size x > 1^1(2. ^^) diverges and is of order (xq — x)"""'^. 
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Construction 


Expected complexity 


C 


= A + B 


hC{x,y) 






C 




AC{x,y) 


= KA{x,y) + KB{x,y) 




C 


= SET>d{B) 


AC{x,y) 


= ''''^>-'-]\,f'''fB{x,y) 


•(l + AS(x,2/)) 


C 


= AolB 


AC{x,y) 


= AAiB{x, y), y) + \A\(^B{x,y),y) ■ AB{x, y) 


C 


= AouB 


AC{x,y) 


^AA{:x,B{x,y)) + \\A\\ 


x,B(x,y)) ■ AB{x,y) 



Figure 13. The expected complexities of Boltzmann samplers specified 
using the samphng rules for the constructions {+,-k,SET>d,°L,°u} (as 
given in Figure [3]) satisfy explicit equations. There exp>_]^(z) = exp(2) 
and, for d> 0, exp>^(z) = X]fe>d z'^/kl- 

Proof. Recall that |C|(a;^j,Q) = x ■ C {x,yo)/C{x,yo), and C is (a — l)-singular if C is a- 
singular. Hence, if a > 1, both functions C(x,yo) and C'{x,yo) converge to positive 
constants as x —^ xq, so that |C|(2;,yQ) also converges to a positive constant. If < a < 1, 
C{x,yo) still converges, but C'{x,yo) diverges, of order (xq — x)°'~^ a.s x ^ xq. Hence 
|C|(3;^yQ) is also of order {xq — a;)"~^. □ 

8.4. Computation rules for the expected complexities of Boltzmann samplers. 

Thanks to Claim [27l the complexity analysis is now reduced to estimating the expected 
complexity AQ"{x,y) when [x,y) is close to a singular point of Q. For this purpose, 
we introduce explicit rules to compute AC{x,y) if C is specified from other classes by a 
decomposition grammar. These rules will be combined with Lemma E5] and Corollarv l26l 
(complexity due to the rejection steps) in order to get a precise asymptotic bound for 
Ag"{x,y). 

We can now formulate the computation rules for the expected complexities. 

Lemma 31 (computation rules for expected complexities). LetC be a class obtained from 
simpler classes A, B by means of one of the constructions {+,*, SET>ci, o^,, o^/}. 

If A and B are equipped with Boltzmann samplers, let VC{x^y) be the Boltzmann sam- 
pler for C obtained from the sampling rules of Figure O Then there are explicit rules, 
as given in Figure \T3[ to compute the expected complexity of TC{x,y) from the expected 
complexities of TA{x,y) andTB{x,y). 

Proof. Disjoint union: TC{x,y) first flips a coin, which (by convention) has unit cost 

in the combinatorial complexity. Then TC{x,y) either calls TA{x,y) or TB{x,y) with 

respective probabilities A{x,y)/C{x,y) and B{x,y)/C{x,y). 

Product: TC{x,y) calls TA{x,y) and then TB{x,y), which yields the formula. 

L- substitution: TC{x, y) calls 7 ^ VA{B{x, y), y) and then replaces each L-atom of 7 by an 

object generated by TB{x,y). Hence, in average, the first step takes time AA[B{x,y),y) 

and the second step takes time \A\(B(x.y).y) ' AB(x,y). 

SET>(i; note that SET>(i(S) is equivalent to A oj^ B, where A := SET>d(Zi), which 
has generating function exp>^(z) := J2k>d '^'^ f^- ^ Boltzmann sampler rA{z,y) simply 
consists in drawing an integer under a conditioned Poisson law Pois>d(z), which is done by 
a simple iterative loop. As the number of iterations is equal to the value that is returned 
(see [TT] for a more detailed discussion) , the expected cost of generation for A is equal to 
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the expected size, i.e. 




Hence, from the formula for A(^ I3){x, y), wc obtain the formula for SET>d. 
U- substitution: the formula for ojj is proved similarly as the one for oj^. 



□ 



Remark 32. When using the computation rules of Figure[T3]in a recursive way, we have 
to be careful to check beforehand that all the expected complexities that are involved 



8.5. Analytic combinatorics of planar graphs. Let C be an a-singular class (see 
Definition [T5|) . A very useful remark to be used all along the analysis of the expected 
complexities is the following: if a > 0, the function C{x,yo) converges when x xq, and 
the limit has to be a positive constant; whereas if a < 0, the function C{x,yo) diverges 
to +00 and is of order {xq — x)" . 

In this section, we review the degrees of singularities of the series of all classes (binary 
trees, dissections, 3-connected, 2-connected, connected, and general planar graphs) that 
are involved in the decomposition of planar graphs. We will use extensively this informa- 
tion to estimate the expected complexities of the Boltzmann samplers in Section [8.61 

Lemma 33 (bicolorcd binary trees). Let TZ = TZ, + TZo be the class of rooted bicolored 
binary trees, which is specified by the system 



Then the classes TZ,, TZo are l/2-singular. The class IC (K,) of rooted (unrooted, resp.) 
asymmetric bicolored binary trees is 1/2-singular (i/ 2- singular, resp.). In addition, these 
two classes have the same singular points as TZ. 

Proof. The classes TZ, and TZo satisfy a decomposition grammar that has a strongly con- 
nected dependency graph. Hence, by a classical theorem of Drmota, Lalley, Woods [Mj, 
the generating functions of these classes have square-root singular type. Notice that, from 
the decomposition grammar Q, the class IC can be expressed as a positive polynomial 
in Zl^ Zjj , TZt, and TZo. Hence IC inherits the singular points and the square-root sin- 
gular type from TZ, , TZo . Finally, the generating function of JC is classically obtained as a 
subtraction (a tree has one more vertices than edges, so subtract the scries counting the 
trees rooted at an edge from the scries counting the trees rooted at a vertex) . The leading 
square-root singular terms cancel out due to the subtraction, leaving a leading singular 
term of degree 3/2. □ 

Lemma 34 (irreducible dissections, from [18]). The class J of rooted irreducible dissec- 
tions is Z/2-singular and has the same singularities as IC. 

Proof. The class is equal to 3 ★ Zl * Zu -kX, which is isomorphic to 3 * Zi^ -k Zjj * /C, so 
J has the same singular points and singularity type as /C. □ 

Lemma 35 (rooted 3-connected planar graphs [2]). The class of edge-rooted 3-connected 
planar graphs is Z 12- singular; and the class of U-derived edge-rooted 3-connected planar 
graphs is 1/2-singular. These classes have the same singular points as IC. 




7^. = Zl*(Z[/ + 7^o)^ tZo (Zu + Tz,f . 
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Proof. The series 63(2, w) has been proved in [24j to have a rational expression in terms 
of the two series R,{z, w) and Ro{z, w) of rooted bicolored binary trees. This property is 
easily shown to be stable by taking derivatives, so the same property holds for the series 
6*3(2:, w). It is proved in [T] that the singular points of G3 are the same as those of TZ, 

and TZo. Hence, the singular expansion of 03(2, w) at any singular point is simply obtained 
from the ones of Rt{z,w) and Ro{z,w); one finds that the square-root terms cancel out, 
leaving a leading singular term of degree 3/2. The study of is similar. First, the rooting 
operator does not change the singular points (as it multiplies a coefficient {n, m) only by 
a factor m), hence, has the same singular points as 7?.,,7?.o, which ensures that the 

singular expansion of 6*3(2, w) can be obtained from those of TZ, and TZo- One finds that 
the leading singular term is this time of the square-root type. □ 

Lemma 36 (networks, from [1]). The classes V, S, V, and TL of networks are 3/2- 
singular, and these classes have the same singular points. 

Lemma 37 (2-connected, connected, and general planar graphs [10]). The classes G2, Gi, 
Q of 2-connected, connected, and general planar graphs are all b/2-singular. In addition, 
the singular points of Q2 are the same as those of networks, and the singular points are 
the same in Qi as in Q . 

8.6. Asymptotic bounds on the expected complexities of Boltzmann samplers. 

This section is dedicated to proving the asymptotic bound KQ"{x, j/o) = 0((.to — x)"^/^). 
For this purpose wc adopt again a bottom-to-top approach, following the scheme of Fig- 
ure (TUl For each class C appearing in this scheme, wc provide an asymptotic bound for the 
expected complexity of the Boltzmann sampler in a neighbourhood of any fixed singular 
point of C In the end wc arrive at the desired estimate of KQ" {x, go). 

8.6.1. Complexity of the Boltzmann samplers for binary trees. 

Lemma 38 (U-derived bicolored binary trees). Let (zo,wo) be a singular point of )C. 
Then, the expected complexity of the Boltzmann sampler for IC — given in Section \^.1.5\ — 
satisfies, 

A^(2, w) = 0^^(20 — 2)""'"/^^ as (2, w) (20, wq). 

Proof. The Boltzmann sampler T!C{z,w) is just obtained by translating a completely 
recursive decomposition grammar. Hence, the generation process consists in building the 
tree node by node following certain branching rules. Accordingly, the cost of generation 
is just equal to the number of nodes of the tree that is finally returned, assuming unit 
cost for building a nodc0. As an unrooted binary tree has two more leaves than nodes, we 
have 

A!C{z,w) < \\JC\\(z,w) < \\K.\\{z,wo), 

where the second inequality results from the monotonicity property of expected sizes 
(Lemma [291). 

Notice that, for t G JC, the number of nodes is not greater than (3|r| -t- 1), where |t| 
is as usual the number of black nodes. Hence the number of nodes is at most 4|r|. As a 
consequence, 

A^(2,?z;) < 4 - 1^1 



We could also use the computation rules for the expected complexities, but here there is the simpler 
argument that the expected complexity is equal to the expected size, as there is no rejection yet. 
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According to Lemma [33l the class JC is 1/2-singular. Hence, by Lemma [501 |K]|(2 i„q) is 
O((zo - z)-^/^) as z ^ zq. So A^(z,w) is also O((zo - z)-i/2). □ 

Lemma 39 (derived bicolored binary trees). Let {zq,wq) be a singular point of IC. Then, 
the expected complexity of the Boltzmann sampler for IC' — given in Section \5.3.1\ — satisfies 

AIC'{z,w) = O ^(zo — z)^"^/^j as {z,w) — > (zo,wo)- 

Proof. The sampler r/C'(z, w) has been obtained from r^(z, w) by applying the procedure 
Uderived— >Lderived to the class IC. It is easily checked that the ratio number of black 
nodes/number of leaves in a bicolored binary tree is bounded from above and from below 
(we have already used the "below" bound in Lemma [38|) . Precisely, 3|r| + 3 > ||r|| 
and |r| < 2||r||/3, from which it is easily checked that a^/u = 2/3 and au/L = 6 
(attained by the tree with 1 black and 3 white nodes). Hence, according to CoroUarv [26l 
A/C'(z,w) < 4A£(z,w), so AK.'{z,w) is O ((zq - z)^^/^). □ 

Lemma 40 (bicolored binary trees). Let {zq,wo) be a singular point of IC. Then, the 
expected complexity of the Boltzmann sampler for IC — given in Section \4.1.6\ — satisfies 

A/C(z, w) = O (1) as (z, w) — > (zq, wq). 

Proof. At each attempt in the generator r/C(z,z/;), the first step is to call r^(z,zi;) to 
generate a certain tree r e £ (it is here convenient to assume that the object is "chosen" 
before the generation starts), with probability 

1 zl-l 



-w 



\r\\. 



K{z,w) \t\\ 

and the probability that the generation succeeds to finish is 2/(||r|| + 1). Hence, the total 
probability of success at each attempt in r/C(z, w) satisfies 



Pace ^ ^ 



1 zl^l ,1^,1 2 



-w 



\\r\\ 



^^^K{z,w) \t\\ \\t\\ + 1 

As each object t ^ K. gives rise to ||r|| objects in K that all have L-size |t| and U-size 
||r|| — 1, we also have 

V- 2 zM^,n,||_i_ 2i^(z,^«) 



Pace = > , -rr, rr-FTW 



A"(z, w) |t|! wK(z, w) 

As IC is 3/2-singular and IC is 1/2-singular, pace converges to the positive constant cq := 
2K{zo,wo)/{woK{zo,wo)) as iz,w) (zo,wo)- 

Now call 2l(z,i(;) the random generator for IC delimited inside the repeat/until loop 
of riC{z,w), and let A^{z,w) be the expected complexity of 2t(z,w). According to 
Lemma [25| AIC{z,w) = A2t(z, i(j)/pacc- In addition, when {z,w) (zo,u'o), Pace con- 
verges to a positive constant, hence it remains to prove that A^{z,w) = 0(1) in order to 
prove the lemma. 

Let T £ IC, and let m := \ \t\\. During a call to 2t(z, w), and knowing (again, in advance) 
that T is under generation, the probability that at least k > 1 nodes of r are built is 
2/(fc + 1), due to the Bernoulli probabilities telescoping each other. Hence, for fc < m — 1, 
the probability pk that the generation aborts when exactly k nodes are generated satisfies 
pi^ = 1^^-^ — = (fc+i)%+2) • addition, the probability that the whole tree is generated 
is 2/m (with a final rejection or not), in which case (m — 1) nodes are built. Measuring 
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the complexity as the number of nodes that are built, we obtain the following expression 
for the expected complexity of ^{z,w) knowing that r is chosen: 

^ — ' m 

k=l 

where Hm := X^fcLi 1/^ is the mth Harmonic number. Define a„i{z) := [w"^]K_{z, w). We 
have 

2 2 

A2l(z,w) < — Viy™a„(z)w™ < — V77„a™(zo)w^. 

K(z,w) ^-^ K[z,w) 

Hence, writing cq := 3/K(z[], wq), we have A2l(z,w) < coj^m -^m'^mi^o)'^™ ^'^^ i^jw) 
close to {zojWq). Using the Drmota-Lalley- Woods theorem (similarly as in Lemma [33)) . it 
is easily shown that the function w i-^ K_{zo, w) has a square-root singularity at w = wq. 
Hence, the transfer theorems of singularity analysis [HI [T3j yield the asymptotic estimate 
o-mizo) cui'^/'^Wq"^ for some constant c > 0, so that am{zo) < c' m~'^/^Wg"'" for some 
constant c' > 0. Hence A2l(z, w) is bounded by the converging scries cq c' ^^^^ Hm mT^I"^ 
for (z, w) close to {zq^wq), which concludes the proof. □ 

8.6.2. Complexity of the Boltzmann samplers for irreducible dissections. 

Lemma 41 (irreducible dissections). Let (zo,wo) be a singular point ofT. Then, the 
expected complexities of the Boltzmann samplers for X and X' — described respectively in 
Section \4.1.T\ and \5.3.2\ — satisfy, as {z,w) (zq^wq): 

AX{z,w) = O (1), 

AX'iz,w) = O [{zo - z)-'/^y 

Proof. As stated in Proposition ! Ill and proved in [18j . the closure-mapping has linear time 
complexity, i.e., there exists a constant A such that the cost of closing any binary tree k is 
at most A • ||k||. Recall that TX{z,w) calls riC{z,w) and closes the binary tree generated. 
Hence 

AXiz, w) < A/C(z, w) + X- ||/C||(,,^) < A/C(z, w) + X- ||/C||(,,„„), 

where the second inequality results from the monotonicity property of expected sizes 
fLemma [29)) . Again we use the fact that, for r G /C, ||r|| < 3|r| -I- 1, so ||r|| < 4|r|. Hence 

AX{z, w) < A/C(2, w)+AX- |/C|(^,„„). 

As the class IC is 3/2-singular. the expected size |/C|(2,miq)) is 0(1) when z Zq. In addi- 
tion, according to Lemma \40\ A/C(z,w) is 0(1) when {z,w) (zq,wo). Hence AX{z,w) 
is 0(1). 

Similarly, for X', we have 

AX'iz, w) < A/C'(z, w) + X- ||/C'||(,,^)) < A/C'(z, w)+4X- |2l */C'| 

As the class K.' is 1/2-singular (and so is Zl *IC'), the expected size \Zl *AC'|(z,u;o) is 
O((zo ~ ^)~^^^) when z — > zq. In addition we have proved in Lemma [39] that A/C'(z, w) is 
O((zo - 2)"^/^). Therefore AX'{z,w) is O((zo - z)"^/^)- □ 

Lemma 42 (rooted irreducible dissections). Let {zo,wo) be a singular point ofX. Then, 
the expected complexities of the Boltzmann samplers for J and J' — described respectively 
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in Section \4-l-T\ and \5.3.2\ — satisfy, as {z,w) (zo,wo)- 

AJ{z,w) = 0(1), 
AJ'{z,w) = O (^{zo - z)-'/^y 

Proof. The sampler TJ^{z, w) is directly obtained from rX(z, w), according to the identity 
J ~ 3-k Zl * Zjj *X, so AJ'{z, w) = AX{z, w), which is 0(1) as (z, w) — > {zq, wq). 

The sampler TJ^'(z, w) is obtained from TX(z, w) and TX' [z, w), according to the iden- 
tity J^' = 3*Zi ★Zc/ *I' + 3*Z,7 Hence, AJ'{z,w) < 1 + AI{z,w) + AI'{z,w). 
According to Lemma HTl Al{z, w) and Al'{z, w) are respectively 0(1) and O{{zo — z)^^^^) 
when {z,w) {zo,wo). Hence AJ'{z,w) is O((zo — 2)""^^^)- D 

Lemma 43 (admissible rooted irreducible dissections). Let {zq^wq) he a singular point of 
X. Then, the expected complexities of the Boltzmann samplers for J7a and J'^ — described 
respectively in Section \4-.l.T\ and \5.3.2\ — satisfy, as {z,w) {z{,,wo): 

AJ^{z,w) = O (1), 

AJi{z,w) = O((zo-^)-'/'). 

Proof. Call TJ{z, w) the sampler that calls VJ{z, w) and checks if the dissection is admis- 
sible. By definition, rj7a(z, w) repeats calling TJ^{z, w) until the dissection generated is in 
jTa- Hence the probability of acceptance Pacc at each attempt is equal to Ja(z, w)/ J{z, w), 
i.e., is equal to Mz{z,w)/ J{z,w) (the isomorphism J7a — A^3 yields Jg_{z,w) = M2,{z,w)). 
Call AJ{z,w) the expected complexity of TJ{z,w). By Lemma 

AMz,w) = AJ{z,w) = _i ' ' AJ{z,w). 

Pace 71/3(0, w) 

We recall from Section 18.51 that the singular points are the same for rooted 3-connected 
planar graphs/maps, for bicolored binary trees, and for irreducible dissections. Hence 
(zq, Wo) is a singular point for M^^z^w). The classes and A^a ~ 2*^/3 are 3/2-singular 
by Lemma [Ml and Lemma [551 respectively. Hence, when {z,w) {zq,wo), the series 
J{z,w) and M3{z,w) are 6(1), even more they converge to positive constants (because 
these functions are rational in terms of bivariate series for binary trees). Hence Pacc also 
converges to a positive constant, so it remains to prove that AJ{z,w) is 0(1). Testing 
admissibility (i.e., the existence of an internal path of length 3 connecting the root-vertex 
to the opposite outer vertex) has clearly linear time complexity. Hence, for some constant 

AJ{z,w) < AJ{z,w)+\-\\J\\(,^^) < AJ{z,w) + X-\\J\\^,,^,-), 
where the second inequality results from the monotonicity of the expected sizes (Lemma l29| . 
Both AJ[z,w) and || J'||(z,t„o) are 0(1) when z z^ (by Lemma [42] and because J is 
3/2-singular, respectively). Hence AJ{z,w) is also 0(1), so AJs,{z,w) is also 0(1). 
The proof for is similar. First, we have 

AJi{z,w) = d^i^^-Aj'{z,w), 
M3'{z,w) 

where AJ'{z,w) is the expected cost of a call to TJ'{z,w) followed by an admissibility 
test. Both series J'{z,w) and M3{z,w) are 1/2-singular, even more, they converge to 
positive constants as {z,w) {zq^wq) (again, because these functions are rational in 
terms of bivariate series of binary trees). Hence, when {z,w) {zq,wq), the quantity 
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J' {z , w) / M3' {z , w) converges to a positive constant. Moreover, according to the linear 
complexity of admissibility testing, we have AJ''{z, w) < KJ\z, -u;) + A • || J''||(z,ujo)- Both 
quantities A (2, ui) and || J''||(2,too) are O((zo — z)^^/^). Hence AJ'jj'(z, u;) isalsoO((zo — 
z)-i/2). ' □ 

8.6.3. Complexity of the Boltzmann samplers for 3-connected maps. 

Lemma 44 (rooted 3-connccted maps). Let (zojWo) be a singular point ofMs. Then the 
expected complexities of the Boltzmann samplers for M.3 and M3' satisfy respectively, as 
{z,w) {zq,Wo): 

kMl{z,w) = 0(1), 
AA^'(2,H = 0({z^-z)-^'^). 

Proof. Recall that TM3,{z,w) {VM3,' {z,w), resp.) calls TJa,{z,w) {rjl{z,w), resp.) and 
returns the primal map of the dissection. The primal-map construction is in fact just a 
reinterpretation of the combinatorial encoding of rooted maps (in particular when deal- 
ing with the half-edge data structure). Hence AMsiz, w) = Aj7a(2, w) and AMs'iz, w) = 
AJ'l{z, w). This concludes the proof, according to the estimates for Aj7a(z, w) and AJ!^{z, w) 
given in Lemma l43l (A proof following the same lines as in Lemma |4T] would also be pos- 
sible.) □ 

8.6.4. Complexity of the Boltzmann samplers for 3-connected planar graphs. 

Lemma 45 (rooted 3-connected planar graphs). Let {zq^wq) be a singular point of Q^. 
Then the expected complexities of the Boltzmann samplers for Q3, Q3' and Q3 satisfy 
respectively, as {z,w) (zo,wo)-' 

Agl{z,w) = 0(1), 

Agl\z,w) = O (^{zo - z)-'/') , 

Agliz,w) = O [{zo - z)-'^^y 

Proof. The sampler TQ3{z,w) (TQ^' {z,w), resp.) is directly obtained from TAi^^z^w) 
(rA^3'(2;, w), resp.) by forgetting the embedding. Hence AG3{z,w) = AM3{z,w) and 
AQ3'{z,w) = AMs'iZjW), which are — by Lemma Hi] — respectively 0(1) and 0((zo — 
z)"i/2) as {z,w) {zo,wq). 

Finally, the sampler TQ3(z,w) is obtained from TQ3'{z,w) by applying the procedure 

Lderived^Uderived to the class Gs- By the Euler relation, ajj/i^ = 3 (given asymp- 
totically by triangulations) and 0^/(7 = 2/3 (given asymptotically by cubic graphs). 

Thus, by Corollary [261 A^(z,w) < 2 • AQ3'{z,w), which ensures that A^(z,w) is 
O{{zo-z)-y''). □ 

8.6.5. Complexity of the Boltzmann samplers for networks. At first we need to introduce 
the following notations. Let C be a class endowed with a Boltzmann sampler TC{x, y) and 
let 7 e C. Then AC^''\x,y) denotes the expected complexity of TC{x,y) conditioned on 
the fact that the object generated is 7. If VC{x,y) uses rejection, i.e., repeats building 
objects and rejecting them until finally an object is accepted, then AC''-' (a;, y) denotes the 
expected complexity of TC{x,y) without counting the last (successful) attempt. 
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Lemma 46 (networks). Let {zo,yo) be a singular point ofD. Then, the expected com- 
plexity of the Boltzmann sampler for T) — described in Section \4-.S\ — satisfies 

KVi^z, yo) = 0(1) as z ^ zq. 

Proof. Trakhtcnbrot's decomposition ensures that a network 7 G I? is a collection of 3- 
connectcd components ki, . . . , (in Q3) that are assembled together in a series-parallel 
backbone (3 (due to the auxihary classes S and V). Moreover, if 7 is produced by the 
Boltzmann sampler r'D{z,yo), then each of the 3-connectcd components Ki results from 
a call to T03{z,w), where w := D{z,yo). 

An important point, which is proved in [T], is that the composition scheme to go from 
rooted 3-connected planar graphs to networks is critical. This means that wq := D{z,yo) 
(change of variable from 3-connected planar graphs to networks) is such that (zq, wq) is a 
singular point of Q^. 

As the series-parallel backbone is built edge by edge, the cost of generating f3 is simply 
||/3|| (the number of edges of f3); and the expected cost of generating ki, for i E [l-.f], is 
A^(''')(z,w). Hence 

r 

(17) AP^-^) {z, yo) = 1 1/3| I + E ^^^"'^ 

i=l 

Claim 47. There exists a constant c such that, for every k Q3, 
Ag3'^'^\z,w) < c\\k\\ as {z,w) ^ {zo,wo). 

Proof of the claim. The Boltzmann sampler rC/3(z, w) is obtained by repeated attempts 
to build binary trees until the tree is successfully generated (no early interruption) and 
gives rise to a 3-connected planar graph (admissibility condition). For k £ IC, call c'-"^ 
the cost of building k (i.e., generate the underlying binary tree and perform the closure). 
Then 

A^(«)(z,H = A^^°j(z,u;) + c(«). 

Notice that AQ/°^{z,w) < AQ3{z,w), which is 0(1) as {z,w) {zo,wo). Moreover, the 
closure-mapping has linear time complexity. Hence there exists a constant c independent 
from K and from z such that AQ^'-'^^z, w) < c\\k\ \ as z zq. A 

The claim ensures that, upon taking c > 1, every 7 G I? satisfies 

r 

Apf^n^, 2/0) <c(||/3||+^ 11^,11) as z^zo. 

1=1 

Since each edge of 7 is represented at most once in /3 U ki U . . . U Kr, we also have 
A£)(T)(z,yo) < c||7||. Hence, when z zq, AP('^)(z,yo) < 3c- (I7I 4- 1) (by the Euler 
relation), which yields 

A2?(z,j/o) < 3c- |Zl*P|(,,^„). 
As the class T> is 3/2-singular (clearly, so is Zl ^T)), the expected size \Zl *2?|(z_j^(,) is 
0(1) when z zq. Hence AI?(z,yo) is 0(1). ' □ 

Lemma 48 (derived networks). Let {zQ,yo) be a singular point ofT>. Then, the expected 
complexity of the Boltzmann sampler for T>' — described in Section \5.5\ — satisfies 

AV'{z, yo) = O ((zo - z)-i/2) as z zq. 
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Proof. Let us fix z e (0,zo). Define X := {AV'{z,yo),AS'{z,yo),AV'{z,yo),An'{z,yo)). 
Our strategy here is to use the computation rules (Figure [T3|) to obtain a recursive equation 
specifying the vector X. By Remarlc [5^ we have to check that the components of X are 
finite. 

Claim 49. For z E (0,zo), the quantities AT)' {z , yo) , AS'{z,yo), AV' {z^yo) , and ATL' [z , yo) 
are finite. 

Proof of the claim. Consider AP'(z, y^) (the verification is similar for AS' {z, yo), AV' {z, yo), 
and ATC'{z,yo)). Let 7 £ V , with (3 the series-parallel backbone and ki, . . . , the 3- 
connected components of 7. Notice that each ki is drawn either by TQ2,{z, w) or TQ3{z, w) 
or rQ3'{z,w), where w ~ D{z,yo). Hence the expected cost of generating is bounded 
by M + c||Ki||, where M Ma,x{AG3{z , w) , AQ3_{z , w) , AQ^' {z , w)) and c||Ki|| represents 
the cost of building Ki using the closure-mapping. As a consequence, 

r 

AV'^'^\z,yo) < ||/3|K^A/ + c||k,|| <C||7||, with C A/ -f c + 1. 

i=l 

Hence 

AI?'(z,2;o)<7577^E IMlfrlyo'"' 
D'{z,yo) |7|! 

which is 0(1) since it converges to the constant CyodyD'{z,yo)/D'{z,yo). A 
Using the computation rules given in Figure (jl3p . the decomposition grammar (N') of 
derived networks — as given in Section [5.51 — is translated to a linear system 

X ^ AX + L, 

where A is a 4 x 4-matrix and L is a 4- vector. Precisely, the components of A are rational 
or exponential expressions in terms of series of networks and their derivatives: all these 
quantities converge as z — > zq because all the classes of networks are 3/2-singular. Hence 
A converges to a matrix Aq as z zq. In addition, A is a substochastic matrix, i.e., a 
matrix with nonncgative coefficients and with sum at most 1 in each row. Indeed, the 
entries in each of the 4 rows of A correspond to probabilities of a Bernoulli switch when 
calling TD'{z, y), TS'{z, y), TP'{z, y), and TH'{z, y), respectively. Hence, the limit matrix 
Aq is also substochastic. It is easily checked that ^0 is indeed strictly substochastic, i.e., 
at least one row has sum < 1 (here, the first and third row add up to 1, whereas the 
second and fourth row add up to < 1). In addition, Ao is irreducible, i.e., the dependency 
graph induced by the nonzero coefficients of Ao is strongly connected. A well known 
result of Markov chain theory ensures that (/ — ^0) is invertible [22]. Hence, (/ — A) is 
invertible for z close to zq, and [I — A)~^ converges to the matrix (/ — Aq)"^. Moreover, 
the components of L are of the form 

L = {a, b,c,d- Ag^'iz, w) + e ■ A^(z, w)j , 

where w = D{z, yo) and {a, b, c, d, e} are expressions involving the series of networks, their 
derivatives, and the quantities {AD, AS, AP, AH}, which have already been shown to be 
bounded as z — > zg. As a consequence, a,b,c,d,e are 0(1) as z — ^ zq. Moreover, it has 
been shown in [T] that the value Wq := D{zo, yo) is such that (zq, wq) is singular for C/3, and 
Wo — w ~ A • (zo — z), with A := D'{zo,yo)- By Lemma l45l AQz'{z,w) and AQz{z, w) are 
O((zo — z)^^/^) as z ^ Zq; hence these quantities are also O((zo — z)~^/^). We conclude 
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that the components of L are 0((zo— z) ^/^), as well as the components oiX ~ {I— A) ^L. 
In particular, AP'(z,yo) (the first component of X) is 0{{zo — z)~^/^). □ 



8.6.6. Complexity of the Boltzmann samplers for 2-connected planar graphs. 

Lemma 50 (rooted 2-connected planar graphs). Let (zo,?;o) be a singular point of Q2- 
Then the expected complexities of the Boltzmann samplers for Q2 and Q2 satisfy respec- 
tively, as z ^ zq: 

A^(z,2;o) = 0(1), 
A^'(z,2/o) - O ((zo - z)-i/2) . 

Proof. Recall that the Boltzmann sampler TQ2{z,yo) is directly obtained from TD{z,yo), 
more precisely from T{l+'D){z, yg). According to Lemma H51 A'D{z, yo) is 0(1) as z ^ zq, 

hence AG2{z,yo) is also 0(1). 

Similarly rQ2'{z, yo) is directly obtained from ri?'(z, j/q)- hence AQ2'{z, y^) — AVlz, j/o), 
which is O((zo — z)~^/'^) as z ^ zp. □ 

Lemma 51 (U-derived 2-connected planar graphs). Let (zg, yo) be a .singular point of Q2. 
Then, the expected complexities of the Boltzmann samplers for Q2 and Q2' — described in 
Section 15.51 — satisfy, as z ^ zo : 



A^(z,2;o) = O (1), 
A^'(z,2;o) = o((zo-z)-i/2). 

Proof. The Boltzmann sampler for ^ is directly obtained from the one for Q2, according 

to the identity 2*^ = * C/2. Hence A^(z,?/o) = A^2(z,?;o), which is 0(1) as 
z — > zq, according to Lemma 1501 Similarly, the Boltzmann sampler for Q2 is directly 

obtained from the ones for the classes Q2 and Q2' ^ according to the identity 2 ★ Q2' = 

* ^' + 2 * Zl ★ 0^2. Hence A^(z, yo) < 1 + A^'(z, yo) + A^(z, yo). When z ^ zq, 

AQ2{z,yo) is 0(1) and AQ2{z,yo) is O((zo — z)~^/^) according to Lemma [5U1 Hence, 
A^'(z,yo)isO((zo-z)-i/2). □ 

Lemma 52 (bi-derived 2-connected planar graphs). Let {zo,yo) be a singular point ofQ2- 
Then, the expected complexities of the Boltzmann samplers for Q2 and Q2 — described in 
Section 15.51 — satisfy, as z —^ zo : 

Ag^{z,yo) = 0(1), 

^G^'{z,yo) = of(zo-z)-i/2 



Proof. Recall that the Boltzmann sampler TQ2{z,yo) is obtained from TQ2{z,yo) by ap- 
plying the procedure Uderived^Lderived to the class C/2. In addition, according to 
the Euler relation, any simple connected planar graph 7 (with I7I the number of vertices 
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and II7II the number of edges) satisfies I7I < II7II + I (trees) and ||7|| < 3|7| — 6 (triangula- 
tions). It is then easily checked that, for the class Q2, ctjj/L = 3 (attained asymptotically 
by triangulations) and a^/^/ = 2 (attained by the link-graph, which has 2 vertices and 1 
edge). Hence, by Corollary [26l KQiiz^yo) < 6A^(z,?/o)- Thus, by Lemma [CTl hQ-iiz^yo) 
is 0(1) as 2: ^ zq. 

The proof for AQ2'{z,yo) is similar, except that the procedure Uderived— >Lderived 
is now applied to the derived class G2' , meaning that the L-sizc is now the number of 
vertices minus 1. We still have au/i^ = 3 (attained asymptotically by triangulations), and 
now a^fu = 1 (attained by the hnk-graph) . Corollarv[26l yields AQ4'(z. vn) < 3 A^'(2;, ?/o)- 
Hence, from Lemma [FTl KQ2 {z,yf^) is 0{{zi^ — z)~^/^) as z — > zq- D 

8.6.7. Complexity of the Boltzmann samplers for connected planar graphs. 

Lemma 53 (derived connected planar graphs). Let (a;o,?/o) a singular point of Qi- 
Then, the expected complexity of the Boltzmann sampler for Qi — described in Section lJ73\ — 
satisfies 

A5i'(x,yo) = 0(1) asx^xo- 



Proof. Recall that the Boltzmann sampler for Qi results from the identity (block decom- 
position. Equation p2)) ) 

gi' = SET {g^OL (Zl*^!')). 

We want to use the computation rules (Figure fO)) to obtain a recursive equation for 
AQi{x, yo). Again, according to Remark[221 we have to check that AQi{x, yo) is finite. 

Claim 54. For < x < xq, the quantity AQi{x,yQ) is finite. 

Proof of the claim. Let 7 G ^1 , with ki, . . . , the 2-connectcd blocks of 7. Wc have 

r 

Ag^^"'\x,yo) = 2II7II +^ Ag/«')(z,2/o), where z = xG^{x,yo). 

1=1 

(The first term stands for the cost of choosing the degrees using a generator for a Poisson 
law; note that the sum of the degrees over all the vertices of 7 is 2| I7I |.) It is easily shown 
that there exists a constant M such that AQ2'^'^'^ {z,ya) < M||k|| for any n € QJ (using 
the fact that such a bound holds for AP^'')(z, j/o) and that rt/2'(z,yo) is obtained from 
ri?(z, yo) via a simple rejection step). Therefore Ag{^'^'> {x, yo) < C| I7I |, with C = 2 + M. 
We conclude that 

-f&yi 

which is 0(1) since it converges to the constant CyodyGi{x,yo)/Gi{x,yQ). A 
The computation rules fFigure [T^]) yield 

^Gi{x,yo) = G2(z, yo) • (Ag2{z,yo) + |52|(z,yo) ' ^Gii^^Vo)) where z = xGi{x,yo), 
so that 

f,r'f ^ Gd{z,yo)Ag-{{z,yo) 

A^i (x, yo = — r— . 

1 - Gd{z,yo) ■ |e?2'|(z,yo) 

Similarly as in the transition from 3-connected planar graphs to networks, we use the 

important point, proved in |20j . that the composition scheme to go from 2-connected 

to connected planar graphs is critical. This means that, when x ^ xq, the quantity 
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z = xG{{x, yo) (which is the change of variable from 2-connected to connected) converges 
to a positive constant zq such that (zo,2/o) is a singular point of Q2- Hence, according 
to Lemma [5^ AQ2{z,yo) is 0(1) as a; ^ xq. Moreover, as the class G2 is 3/2-singular, 
the scries G'2'(^,J/o) Stiid the expected size |52'|(z,i(o) converge to positive constants that 
are denoted respectively G2 {zQ^y^) and \Q2\{za,yo)- have shown that the numerator 
of AQi {x^yo) is 0(1) and that the denominator converges as a; ^ xq. To prove that 
KQ{{x, yo) is 0(1), it remains to check that the denominator does not converge to 0, i.e., 
to prove that Gdizo.yo) ■ \Q2\(zo,ya) 7^ 1- 

To show this, we use the simple trick that the expected complexity and expected size 
of Boltzmann samplers satisfy similar computation rules. Indeed, from Equation (|12|1. it 
is easy to derive the equation 

\Q(\(x,yo) = G2{z,yo) ■ \Q2\(z,yo) ' {\Ql\(x.ya) + l) whcrc z = xG{{x,ya), 

either using the formula |C|(j;.y) = dxC{x,y)/G{x,y), or simply by interpreting what 
happens during a call to TQi{x,y) (an average of G2{z,yo) blocks are attached at the 
root- vertex, each block has average size |52'|(z,i;o) ^^'^ carries a connected component of 
average size {\Gi\{x,yo) + 1) ea,ch non-root vertex). Hence 

|^,| G2iz,yo)-\g2\{z,yo) 

'^'''"^"^~l-G2'(.,2/o)-|e2'|(,,o)- 

Notice that this is the same expression as Agi'{x,yo), except for |^2'|(z,j,o) replacing 
(-2, yo) in the numerator. The important point is that we already know that |Si |(a;,yo) 
converges as x ^ xo, since the class g{ is 3/2-singular (see Lemma [57| . Hence G2{zo, yo) ■ 
\g2'\{zo,yo) to be different from 1 (more precisely, it is strictly less than 1), which 
concludes the proof. □ 

Lemma 55 (bi-derived connected planar graphs). Let {xo,yo) be a singular point of 
Qi. Then, the expected complexity of the Boltzmann sampler for Q\ — described in Sec- 
tion \5.5.1\ — satisfies 

Agi"{x,yo) = O ((a-o - x)-i/2) as x ^ xq. 



Proof. The proof for AQi" {x,yo) is easier than for AQi (x^ya). Recall that TQi" {x,yo) is 
obtained from the identity 

Si" = (e/ + Zi * ^i") * g^' OL [Zl * Q{) * Q{. 

At first one easily checks (using similar arguments as in Claim [54)) that AQi{x, yo) is finite. 
Using the computation rules given in Figurc [T51 we obtain, writing as usual z = xG{{x, yo), 

\r"( \ 11 G{{x,yo) . , xG{'{x,yo) „ , 

Gi'{x,yo)+xGf[x,yo) Gi'{x,yo)+xGf{x,yo) 
+A(?2"(2, yo) + \G2\(z.y„) ■ Ag{{x, yo) + AQ{{x, yo). 

Hence 

^5iix,yo) = aix,yo) ■ {l + b{x,yo) ■ Ag^{x,yo) + Ag^'{z,yo) + \g^'\(^,,y„y AQ^ix^yo)), 
where 

Gi{x,yo) + xGi'{x,yo) - 2Gi{x,yo) + xGi'{x,yo) 

a[x,yo) = 7777 r , 02^,2/0) = -7777 T"^ — F^tt r- 

G]!ix,yo) G{[x,yo) + xGf{x,yo) 



UNIFORM RANDOM SAMPLING OF PLANAR GRAPHS IN LINEAR TIME 



53 



As the classes g{ and Gi are respectively 3/2-singular and 1/2-singular, the series a{x, yo) 
and b{x,yo) converge when x xq. As Q2' is 1/2-singular, \G2"\{z,yo) C>{{zo — z)~^/^) 
when z zq. Moreover, according to Lemma \52\ AQ2" {z,yo) is O{{zo — z)^^!"^). Next 
we use the fact that the change of variable from 2-connected to connected is critical. 
Precisely, as proved in [1], when a; — ^ and when z and x are related by 2; = xG{{x, yo), 
we have zq ~ z ^ X ■ [xq — x), with A := limdz/dx = xqGi" {xo,yo) + Gi {xo,yo). Hence, 
\G2'\(z,yo)) ^^'^ AC/2"(z,yo) are O{{xo — x)~^^^). In addition, we have proved in Lemma [53l 
that Agi{x,yo) is 0(1). We conclude that Agi{x,yo) is O{{xo - □ 

Lemma 56 (connected planar graphs). Let {xo,yo) be a singular point oj Q\. Then, the 
expected complexity of the Boltzmann sampler for Qi — described in Section \4.!^ — satisfies 

At/i(a;, j/o) = O (1) as a; — > xq. 

Proof. As described in Section [43l the sampler TQi{x,y) computes 7 ^ TQi{x,y) and 
keeps 7 with probability 1/(|7| + 1). Hence the probability of success at each attempt is 

^ 1 1 a^'^' II7II ^ 1 a^'^' II7II 

^^^^ G{{x,yo)^^,\l\ + l\lf' G{{x,yo)^^,{\l\ + l)\^' ' 

Recall that for any class C, ,„ identifies to Cn+i.m- Hence 

_ 1 _ Gi{x,yQ) 

P-^^- G{{x,y,) ^ |7|! ~ xG{{x,y,y 

In addition, by Lemma 1251 A^i(.t,?;o) = ^Qi {^tVo) I Pace- As the classes Qi and Qi are 
respectively 5/2-singular and 3/2-singular, both series Gi(a;,yo) and Gi{x,yo) converge 
to positive constants when a: —> xq. Hence pacc converges to a positive constant as well. 
In addition. At// (a;,?/o) is 0(1) by Lemma [53l Hence AC/i(x,yo) is also 0(1). □ 

8.6.8. Complexity of the Boltzmann samplers for planar graphs. 

Lemma 57 (planar graphs). Let (.to,?/o) be a singular point ofQ. Then, the expected com- 
plexities of the Boltzmann samplers for Q , Q' and Q" — described in Section \4-.4\ and \5.5.^ — 
satisfy, as x Xq: 

Kg{x,yo) = 0(1), 
Ag'{x,yo) = 0(1), 
Kg"{x,yo) = O ((xo-x)-i/2). 

Proof. Recall that TQ(x,y) is obtained from TQi{x,y) using the identity 

9 = SET(gi), 

hence KQ{x,yQ) = Gi{x,yt)) ■ KQi(x,yt)). When a; — > xq, Gi{x,yo) converges (because Gi 
is 5/2-singular) and AQi{x,yo) is 0(1) (by Lemma [55)1 . Hence AQ{x,yo) is 0(1). 
Then, TG'{x,y) is obtained from TQ{{x,y) and TQ(x,y) using the identity 

G' ^Gi^g. 

Hence A^;'(a;,yo) = Agi'{x,yo) + Agix,yo). When x xq, A^i'(x,yo) is 0(1) (by 
Lemma [531) and Ag{x,yo) is 0(1), as proved above. Hence Ag'{x,yo) is 0(1). 

Finally, rQ"{x,y) is obtained from rQ{'{x,y), TQ{{x,y), TQ'{x,y), and TQ{x,y) using 
the identity 
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Hence 

Ag"(.T, yo) = 1 + ^ (A5i"(x, yo) + m{x, yo)) + ^ (A5i'(x, y^) + Ag'(x, yo)) , 
a + a + 

where a = {x,yo)G{x,yo) and 6 = G{(x,yt))G' {x,yQ). Thus 

A5"(a-,yo) < l + Agl"(.T,yo)+Ag(.T,yo)+Aal'(2;,2/o)+Ae'(x,^/o)• 
When a; ^ xo, A^;i"(x,?/o) is 0((a-o - x)'^''^) (by Lemma [53, A^;i'(x,?;o) is 0(1) (by 
Lemma [55]) ■ and AQ'{x,yo) and AQ{x,yo) are 0(1), as proved above. Hence AQ"{x,yo) is 
O{{xo — which concludes the proof. □ 

This concludes the proof of the expected complexities of our random samplers. (Recall 
that, thanks to Claim [27l the proof has been reduced to proving the asymptotic estimate 
Ae"(x,2/o)=0((xo-x)-i/2).) 
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